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Study guide 


There are five sections to this chapter. They are 
intended to be studied consecutively in five study 
sessions. 


Section 1 requires the use of an audio CD player, 
and Section 5 requires the use of the computer 
together with Computer Book B. 


The pattern of study for each session might be as 
follows. 


Study session 1: Section 1. 
Study session 2: Section 2. 
Study session 3: Section 3. 
Study session 4: Section 4. 
Study session 5: Section 5. 


Each session will require two to three hours to 
study, the longest being the second. 


In addition, there is an optional video band 
associated with this chapter, in which an 
illustration is given of the possible uses of certain 
linear transformations in design. The best time to 


watch this band is after you have studied Section 4. 


Before studying this chapter, you should be 
familiar with the following topics: 


© the concept of a fixed point of a function; 


© linear transformations, in particular, rotations, 
reflections, scalings and shears; 


© the matrix representation of a linear 
transformation; 


© the determinant test for the invertibility of a 
2 x 2 matrix; 


© matrix algebra including matrix multiplication 
and the formula for the inverse of a 2 x 2 
matrix, involving its determinant; 


© the long-term behaviour of the geometric 
sequence r”, and the evaluation of limits of 
ratios involving geometric sequences. 


The optional Video Band B(iii) Algebra 
workout — Eigenvectors could be viewed at 
any stage during your study of this chapter. 


Introduction 


In this chapter two concepts, iteration of functions from Chapter B1 and 
linear transformations from Chapter B2, are drawn together into the study 
of the iteration of linear transformations. 


The main theme is the further study of linear transformations, but affine 
transformations, which were introduced in Chapter B2, are also 
investigated briefly. 


Any function f with rule x +> Ax, where A is a 2 x 2 matrix and x isa 
vector, represents a linear transformation. You will see that a geometric 
understanding of such functions provides a key to understanding the 
long-term behaviour of their iterations. In particular, it turns out to be 
important to identify certain sets of points that remain unchanged when a 
linear transformation is applied. 


This chapter begins by studying the geometry of certain types of linear 
transformations, and looking at the points which remain fixed under these 
transformations. We also look at lines which remain unchanged under 
these transformations; on such lines, the individual points might not be 
fixed by the transformation, but the image of the entire line is the same as 
the original line. We call such a line an invariant line. Invariant lines 
through the origin are of particular interest. In Section 1, we look at these 
invariant lines from a geometric perspective for certain types of linear 
transformation, but we soon find that we need to investigate them 
algebraically. Section 2 takes an algebraic approach which allows us to find 
such invariant lines for many linear transformations represented by 2 x 2 
matrices. 


In Section 3, you will see how to write certain 2 x 2 matrices A in the form 
PDP '," where D is a diagonal 2 x 2 matrix. This process, called 
diagonalisation, gives a neat and efficient way of finding powers of a matrix. 


In Section 4, we use the results from Sections 2 and 3 to explore the 
iteration of linear transformations. It turns out that the long-term 
behaviour of sequences of points generated by iteration of linear 
transformations can often be described by using properties of the 
corresponding matrix. 


In Section 5, the computer is used to investigate patterns in the iteration 
of linear and affine transformations. A description is given of how iteration 
can lead to pictures of natural objects such as ferns. 


1 Fixed points and invariant lines 


Recall that a linear 
transformation maps the 
origin to itself and maps 
straight lines to straight lines. 


Recall from Chapter B1, 
Section 1, that a fixed point 
of a real function f is a real 
number a such that f(a) = a. 


To study this section you will need an audio CD player and CDA5493. 


When a typical linear transformation is applied to the plane, most points 
are moved to a new position. As a consequence, the geometrical properties 
of a figure (shape, angles, distances, etc.) are altered when the linear 
transformation is applied. However, certain things are not altered by the 
application of a linear transformation — for example, the origin is 
unchanged, and some straight lines may also remain unchanged. In this 
section we investigate points and lines which are unchanged under certain 
basic linear transformations. 


In the audio band, we use our geometric understanding of rotations, 
reflections, scalings and shears to determine the points and lines which 
remain unchanged when these types of linear transformation are applied. 


A point which remains unchanged under a linear transformation is called a 
fixed point of that transformation. The first half of the audio band 
considers fixed points of particular linear transformations. The second half 
of the audio band considers lines which remain unchanged under the same 
linear transformations. These lines are called invariant lines. 


Now listen to CDA5493 (Tracks 11-18), band 4, ‘Fixed points and 
invariant lines’. 
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Fixed Foints 


A fixed point of a linear transformation is a point which is equal to ite image under the 
linear transformation. 
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Fixed points of the rotation rg: 
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Fixed points of the reflection og: 
All points on the axis of reflection 
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Scalings 
= scaling yh 


with 


factors 
3 and -1 


ee 
e 


Fixed points of the scaling with factors 4 and b: 


Factors | a=1, b=1 | a=1 beet | azi,b=1 | ac, bel 


Fixed points | every point | x-axis | y-axis | O 


All points 
move 
except the 
origin. 


Frame 6 


Solutions are given on page 52. 
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Invariant lines 


An invariant line of a linear transformation is a line that is equal to its image line. 
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Invariant lines through the origin for the reflection 4: 
Two invariant lines — the axis of reflection and the line perpendicular to the axis of reflection. 
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Activity 1.2 Invariant lines through the origin 


Describe the invariant lines through the origin (if there are any) of each of 
the following linear transformations. 


(a) T3n/2 (b) dz/2 (c) scaling with factors 1 and —1 
(d) y-shear with factor 4 


Solutions are given on page 52. 


The fixed points and invariant lines through the origin for each of the four 
types of linear transformation considered in the audio tape are summarised 
below. 


Fixed points 


© Most rotations about the origin move every point except the origin. 
The only exceptions are rotations through an integer multiple of 27, in 
which case every point is a fixed point. 


© A reflection in a line through the origin has exactly one line of fixed 
points: the axis of reflection. 


© Most scalings have no fixed points except the origin. The exceptions 
occur when one (or both) of the factors is one. 


© Most shears have exactly one line of fixed points: for an x-shear this 
line is the x-axis, for a y-shear this line is the y-axis. The exceptions 
are shears with factor zero, for which every point is a fixed point. 


Invariant lines through the origin There may also be other 
invariant lines not passing 


© Most rotations about the origin have no invariant lines. The a 
through the origin. 


exceptions are rotations through an integer multiple of 7, in which 
case every line through the origin is an invariant line. 


© For any reflection in a line through the origin, the axis of reflection is 
an invariant line consisting of fixed points. The line (through the 
origin) perpendicular to the axis of reflection is also an invariant line. 


© A scaling with factors a and 6} has at least two invariant lines: the Remember that the scale 
coordinate axes. However, if a = 6 then every line through the origin factors a and 6 are both 
is an invariant line. non-zero. 


© Most shears have just one invariant line through the origin: for an 
a-shear this line is the x-axis and for a y-shear this line is the y-axis. 
However, if the shear has factor 0, then every line through the origin is 
an invariant line. 


We have used our geometric understanding of these four types of linear 
transformation to determine their fixed points and invariant lines through 
the origin. If we know the matrix which represents one of these linear 
transformations, then we can check our results using algebra. 


For example, reflection in the z-axis is represented by the matrix See Chapter B2, 
1 0 : : ; Subsection 1.3. 

Qo = ( i a }. From our geometric understanding of reflections, we Ss 

know that every point on the z-axis (the axis of reflection) is a fixed point. 

To check this algebraically, we find the image of an arbitrary point on the 


wh 


Indeed, any line 
perpendicular to the z-axis is 
an invariant line. 
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x-axis, say (c,0). We represent the point (c,0) by the vector x = ca and 
show that Ax = x, where A = Qp. In fact, 


1 0 ey Fee). Fe 
0 -1 0/ \0+0/ \O/’ 
so (c,0) is a fixed point, as expected. Since (c,0) is an arbitrary point on 


the x-axis, we have shown that every point on that axis is a fixed point. 


We can define a fixed point of a linear transformation algebraically as 
follows. 


Algebraic definition of a fixed point 


A fixed point of a linear transformation represented by the 
matrix A is a point, represented by the vector x, such that Ax = x. 


Similarly, from our geometric understanding of reflections, we know that 
for the above reflection, the y-axis is an invariant line; that is, that the 
image of each point on the y-axis lies on the y-axis. We now use algebra to 
check this fact. Let (0,c) be an arbitrary point on the y-axis. The image of 
the point (0,c) is given by 


CevOne 


The point (0,—c) lies on the y-axis. Since we started with an arbitrary 
point on the y-axis, we have shown that every point on the y-axis has its 
image on the y-axis. Also, as c varies, this image ranges over the whole 
y-axis. Thus the y-axis is an invariant line for this reflection. 


Activity 1.3 Checking a fixed point and an invariant line 


Use the matrix that represents a scaling with factors 1 and 3 to show 
algebraically that: 


(a) the point (2,0) is a fixed point of this scaling; 
(b) the y-axis is an invariant line of this scaling. 


Solutions are given on page 52. 


Given a line and the matrix representing a linear transformation, we can 
check whether that line is an invariant line of the linear transformation as 
in Activity 1.3. In Section 2, we address the problem of determining such 
invariant lines. 


Summary of Section 1 


This section has introduced: 


© the ideas of a fixed point and an invariant line of a linear 
transformation; 


© the determination, by geometric visualisation, of the fixed points and 
invariant lines through the origin for rotations, reflections, scalings and 
shears; 


© a matrix method for checking such determinations. 
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Exercises for Section 1 


Exercise 1.1 


Describe the fixed points and invariant lines through the origin of each of 
the following linear transformations. 


(a) Tr (b) dna (c) scaling with factors 2 and —3 
(d) a-shear with factor 1 


Exercise 1.2 


This exercise concerns the x-shear with factor 1. 


(a) Use the matrix representing this shear to check that every point on the 
x-axis is a fixed point of this shear. 


(b) Show that the line y = x is not an invariant line of this shear. 
(Hint: Find a point on the line y = x« whose image does not lie on the 
line y= x.) 
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In the previous section, we used our knowledge of the geometric effect of 
certain types of linear transformation to investigate their fixed points and 
invariant lines through the origin. Given the matrix which represents one 
of these linear transformations, you saw how to check such geometric 
results using algebra. However, for a linear transformation represented by 
a general 2 x 2 matrix, you might not have a clear understanding of the 
geometric effect of the linear transformation on the plane. For such a 
linear transformation, the algebraic definition of a fixed point, given in 
Section 1, can be used to identify any fixed points of the linear 
transformation. Identifying the invariant lines through the origin is not so 
straightforward, and in this section we develop an algebraic method to do 
this. Although identifying the fixed points of a linear transformation does 
give us some information about its geometric effect, identifying the 
invariant lines through the origin gives us far more information. 


Throughout the remainder of this chapter, we often consider points in the 
plane R’ and use the vector representation of these points in algebraic 
manipulations. When we talk about transformations of the plane in terms 
of matrices, the point (x,y) will be represented by its position 


vector & . Generally, the word ‘point’ is used when discussing the plane, 


and the word ‘vector’ is used when discussing the related algebra. But to 
avoid cumbersome expressions when both contexts are relevant, the words 
point and vector will often be used interchangeably. 


2.1 Finding eigenvalues and eigenlines 


As a first step towards an algebraic method for finding the invariant lines 
through the origin for a linear transformation, we need to look more 
closely at the images of points on such invariant lines. 


Activity 2.1 Finding invariant lines 


Consider the linear transformation 


f(x) = Ax, where A= € a. 


(a) By finding the image of the arbitrary vector & on the line y = —2, 
show that this line is an invariant line of the linear transformation f. 


(b) Find the image under f of each of the following vectors, which lie on 


“OQ 6 


Express each image as a scalar multiple of the vector whose image it is. 
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(c) By finding the image of the arbitrary vector ( 5° on the line y = 3x, 
2 
show that this line is an invariant line of the linear transformation f. 


(d) Find the image under f of each of the following vectors, which lie on 


the line y = 3a. 
G) @) ©) 
3 3 —3 


Express each image as a scalar multiple of the vector whose image it is. 


Solutions are given on page 52. 


In Activity 2.1(a) and (c), you saw that the lines y = —a and y = x are 
invariant lines of the linear transformation f/f. 


In Activity 2.1(b), you saw that the image of each given vector is —1 times 
the given vector. For example, 


r(3)=(3)=-1(). 


In general, the image of the vector (_£) is (~£) =-l e as shown 
in Figure 2.1. 


Also, in Activity 2.1(d), you saw that the image of each given vector is 4 
times the given vector. For example, 


#(a)=(6)=4(4). 


In general, the image of the vector ( 3 


role 
SY 
a 
n 
aS 
om 
a 0 
Wo 
II 
oe 
aS 
tole 
a 0 
SF 
ie) 
na 
n 
> 
° 
= 
5 
or 
= 


Figure 2.1. 


. . . . . __ — ai 
Figure 2.1 The effect of f on points on the invariant lines y = —a@ and y = 5a 


As indicated in Figure 2.1, each point on an invariant line through the 


origin is mapped to a point which is a fixed scalar multiple of itself. In A scalar multiple of a point 
fact, as is shown below, this property holds for any invariant line through (a, y) has the form (cz, cy), 
the origin of any linear transformation. where ce R. 
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Suppose that the linear transformation represented by the matrix A has 
an invariant line through the origin and that (x, y) is a point on this line, 
represented by the vector x. Then, for some k € R, 


Ax = hx, (2.1) 


Now what will be the effect of the transformation on another point on this 
invariant line? If a point is on the same line through the origin as the 
point (x,y), then it is a scalar multiple of (x, y), say (cx,cy), where cE R. 
So what happens to the point (cx, cy), represented by the vector cx, when 
we apply the transformation represented by A? We obtain 


A(cx) =c(Ax), | by the rule for scalar multiplication of matrices, 
=c(kx), by equation (2.1), 
= k(ex). 


So every point (cx, cy) on this invariant line is also scalar multiplied by the 
number k when the transformation is applied. Any line through the origin 
with this scaling property is called an eigenline and the corresponding 
constant scale factor & is called an eigenvalue of the matrix A. Thus we 
have shown that every invariant line through the origin of a linear 
transformation is an eigenline of the matrix that represents the 
transformation. Any vector representing a non-zero point on the eigenline 
is called an eigenvector. 


1 2 
3 2 
A is y = —a, with corresponding eigenvalue —1 and another eigenline of A 
isy= Sy, with corresponding eigenvalue 4. Each eigenline has many 
eigenvectors; in fact, any non-zero point on an eigenline corresponds to an 
eigenvector. Two eigenvectors for the eigenline y = —x are the vectors 


1 —2 : : : 
(_1) and ( 9 ). Two eigenvectors for the eigenline y = aa are the 


For example, for the matrix A = ( ) in Activity 2.1, one eigenline of 


vectors (3) and = ip These eigenvectors are illustrated in Figure 2.2. 


Points on 


Points on : 
this this 
elgenline eigenline 
are are 
scalar scalar 
multiplied multiplied 
by the by the 


eigenvalue 4 


eigenvalue —1 


Figure 2.2 Eigenvectors 
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Here are the formal definitions of the terms introduced above. 


Definitions 


Let A be a 2 x 2 matrix. If x is a non-zero vector representing a 
point (x,y) for which there is a real number k such that Ax = kx, 
then 


© kis called an eigenvalue of A and x an eigenvector of A; 


© the line through the origin on which an eigenvector lies is called 
an eigenline; 


© Ax = kx is called the eigenvector equation. 


A number of remarks about these definitions are in order. 


ils 


If k is an eigenvalue of a matrix A that represents a linear 
transformation f, then k is often called ‘an eigenvalue of f’. Also, the 
phrases ‘eigenvector of f’ and ‘eigenline of f’ are in common use. 

Not all matrices have eigenvalues. See Activity 2.8, for example. 

The condition that the eigenvector x is non-zero is important. If x is 
the zero vector, then the equation Ax = kx holds for any number k 
and any 2 x 2 matrix A. We are interested in picking out the non-zero 
vectors x such that Ax is a scalar multiple of x, and the values of k 
that appear as ‘scale factors’. 

For each eigenline of A there is a corresponding value of k. For each 
eigenvalue k, the eigenvectors that satisfy the eigenvector equation are 
said to correspond to that value. 

It follows from the eigenvector equation, Ax = kx, that the linear 
transformation f represented by A has a certain scaling effect. 

Figure 2.3 shows this effect on eigenvectors u, and uy, along the two 


eigenlines, €; and 2, of A. In this figure, the corresponding eigenvalues 


k, and ky are taken as positive and negative, respectively. 


yA yA 


Figure 2.3 Scaling along eigenlines 


In general, if (x, y) lies on an eigenline with eigenvalue k, then the 
image of (x,y) is |k| times further from the origin along the eigenline 
than (x,y) is. If k > 0, then the image lies on the same side of the 


origin as (a, y). If k < 0, then the image lies on the opposite side of the 


origin from (2, y). 


The word ‘eigen’ is the 
German adjective meaning 
‘characteristic’ or ‘own’, so an 
eigenvalue of a matrix is a 
number that is characteristic 
of the matrix. 

Many texts use the Greek 
letters A (lambda) and yu 
(mu) to represent eigenvalues. 


Note that each eigenline 
includes the origin, but that 


the vector : is not an 


eigenvector. 


iwi 


This linear transformation is 
not invertible because it is 
not one-one. Hence the 
matrix representing it is not 
invertible. 


The left-hand matrix is 
A — KI, where I is the 2 x 2 
0 


: ‘ ‘ 1 
identity matrix 0] 


The solution of matrix 
equations such as 

equation (2.3) was discussed 
in MST121 Chapter B2, 
Subsection 5.3. 
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6. You will see later that a matrix can have eigenvalue 0. The linear 
transformation represented by such a matrix maps all points on the 
eigenline corresponding to eigenvalue 0 to (0,0). In this case, therefore, 
the eigenline is not an invariant line (since its image is just a point). 
However, if A is an invertible matrix that represents an invertible 
linear transformation f, then the eigenlines of A and the invariant 
lines through the origin of f are identical. 

7. Equations of the form Ax = kx arise in a large variety of mechanical 
and electrical systems in engineering. For example, the resonant 
frequencies of bridges, aeroplanes, cars and other mechanical structures 
are calculated by solving an eigenvector equation. The eigenvector 
equation also arises in statistical, numerical and other non-engineering 
situations, such as the modelling of the growth of populations, as you 
will see in Section 4. 


Finding eigenvalues 


If we are given the matrix which represents a linear transformation, how 
can we determine its eigenlines and eigenvalues (if there are any)? We 
start by looking at an example, before dealing with the general case. 
Suppose we want to find the eigenlines and eigenvalues of the matrix 
ne € ) 
3 2 
invariant line through the origin. Then Ax = kx for some k € R; that is, 


(3 2) (5) =*(): 


This matrix equation corresponds to the two simultaneous equations 


. Let x be a non-zero vector representing a point (2, y) on an 


x+y = ka, 
3a + 2y = ky; 
that is, 


(1—k)x+2y=0, 
32 + (2—k)y =0. (2.2) 


These are two equations in three unknowns zx, y and k. We could begin by 
trying to solve for x and y in terms of k, but in fact it is easier to begin by 
finding k. 


Putting equations (2.2) into matrix form, we obtain 
1-k 2 ee ie 
3 2=h)\y) \0)- 
—k 
3 2-k 


(2.3) 


: 1 — : 
Suppose now that the matrix ( ) is invertible. Then we 


obtain the solution 


GCs" ao 


However, we are looking for solutions of equation (2.3) in which ) is 


—k 2 


3 2 2 must be non-invertible. 


non-zero. So the matrix ( 
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2 


; 1—k 
Hence the determinant of ( 3 o_k 


) must be zero; that is, 
(1—k)(2—k) -6=0. 


Simplifying this equation, we obtain k? — 3k — 4 = 0, which can be solved 
by factorising or by using the quadratic formula to give the eigenvalues 
k =4 and k = —1 (which agree with Activity 2.1). 


Once we have found the eigenvalues of a matrix, it is possible to use these 
values to find the equations of the corresponding eigenlines. We shall 
return to this task after we have discussed finding eigenvalues in the 
general case. 


By using the method just described to find the eigenvalues for a general 
2 x 2 matrix, we shall find a formula which will allow us to calculate 


eigenvalues quickly. Suppose that we want to find the eigenvalues of the 
matrix A = ¢ . 
c d 


(x,y) on an eigenline. Then Ax = kx for some k € R; that is, 


(2 a) (G) ="): 


This matrix equation corresponds to the two simultaneous equations 


iF Let x be a non-zero vector representing a point 


ax + by = ka, 

ca + dy = ky; 
that is, 

(a—k)x+ by =0, 

cx + (d—k)y =0. (2.4) 


Putting equations (2.4) into matrix form, we obtain 


eo ata) GG) = (a) 
c d—-k y 0)° 
Since we want the vector x to be non-zero, the matrix 
(“ —-k b ) 
c d—-—k 
must be non-invertible; that is, it must have zero determinant. So 
(a—k)(d—k) — bc = 0. 
Rearranging this equation, we obtain 
k? —(a+d)k+ad—bc=0. (2.5) 


This quadratic equation is called the characteristic equation of A. 


Thus, for a 2 x 2 matrix A = t a the eigenvalues of A are the 


solutions (if any) of the characteristic equation of A. This equation enables 
us to quickly calculate the eigenvalues of a 2 x 2 matrix — if there are any. 
If the characteristic equation has no real solutions, then A has no 
eigenvalues. 


The determinant test for 
non-invertibility of a matrix, 
based on showing that the 
determinant of the matrix 

is 0, was stated in MST121 
Chapter B2, Subsection 5.2. 
Recall that the determinant 


of the matrix : e is 


d 
ad — be. 


It may help you to remember 
this equation if you notice 
that the coefficient of k is 
minus the sum of the diagonal 
entries of A, called the trace 
of A, and that the constant 
term is the determinant of A. 
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Here the trace is 6 and the 
determinant is 5, so you can 
write down 


k?-6k+5=0 
directly. 
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A handy way to check that you have correctly calculated the eigenvalues of 
A= (: 1) is to check that the sum of the eigenvalues is equal to the 
trace a+d, the sum of the diagonal elements of the matrix A. For 

oS has eigenval —1 
3 9) has eigenvalues 


and 4. Here, the sum of the eigenvalues is —1 + 4 = 3, and the trace is 
a+d=1+2=83 as expected. 


example, you have seen that the matrix A = ( 


When finding the eigenvalues of a particular matrix, there is no need to 
repeat the argument deriving the characteristic equation given above. Just 
write down the characteristic equation, and solve it by factorising or by 
using the quadratic formula. 


Example 2.1 Finding eigenvalues 


Find the eigenvalues (if they exist) of each of the following matrices. 


@ (Gi) (7 3) 


Solution 


3.4 
k? —(2+4)k+8-3=0; 
that is, k? —6k+5=0. This equation can be factorised: 
(k —5)(k-1) =0, 


(a) We have (: a = G , ), so the characteristic equation is 


so the eigenvalues are 1 and 5. (A quick check shows that the sum of 
the eigenvalues, 1 + 5 = 6, is equal to the trace, 2+ 4 = 6.) 


— 
= 


d —1 8 
k?-k+1=0. 


We have (< i) = (= : ), so the characteristic equation is 


If we try to solve this quadratic equation using the formula, then we 
obtain 


k= 3 (14 (-1)? = 4x1). 


Now (—1)? —4 x 1 = —3 is negative, so this quadratic equation has no 
real solutions. Hence this matrix has no eigenvalues. 


Now it is your turn to find some eigenvalues. 


Activity 2.2 Finding eigenvalues 


For each of the following matrices, write down the characteristic equation 
and use it to find the eigenvalues (if they exist). 


@ (22) wm (2 3) © (05 33) 


Solutions are given on page 52. 
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Finding eigenlines 


Now that we know how to find the eigenvalues for a 2 x 2 matrix (when 
they exist), how can we use these eigenvalues to find the corresponding 
eigenlines? We start by returning to our initial example, involving the 


matrix A = e ‘ 


3 9 ). We found that this matrix has eigenvalues k = —1 
and k = 4. 


To find the equation of the eigenline corresponding to the eigenvalue 
k = —1, we substitute k = —1 into the eigenvector equation Ax = kx. 
This gives 


(3 2) (G)=4G), 


This matrix equation corresponds to the following two equations 


r+ 2y=—2, 
3x4 + 2y = —-y; 
that is, 
2x2 + 2y =0, 
3x + 3y = 0. 


These equations both reduce to x + y = 0. Thus the eigenline 
corresponding to the eigenvalue k = —1 has equation y = —2x. 


We repeat this procedure with the eigenvalue k = 4 to find the other 
eigenline. We have 


(2) G)-4G). 
This matrix equation corresponds to the following two equations 
z+ 2y = 42, 
3x + 2y = Ay; 
that is, 
—3x + 2y =0, 
3x — 2y =0. 


These equations both reduce to 3x — 2y = 0. Thus the eigenline 
corresponding to the eigenvalue k = 4 has equation y = Sy. 


Notice in this example that, of the pair of equations obtained from the 
eigenvector equation, one is redundant each time. This always happens 
(provided that both equations involve both xz and y). This fact provides a 
useful practical check on your working: having found the eigenline using 
one of the equations, you can check that the other equation is satisfied too. 
If not, you have made a mistake! 


Once we know the equations of the eigenlines, we can choose eigenvectors 
for each eigenline. An eigenvector for the eigenline y = —x is any non-zero 


vector u = (*) for which the corresponding point (u,v) lies on the 


2 : : F F 1 : 
eigenline. In this case, one eigenvector is & ), Two other eigenvectors are 


(<2) (Vs) 


Jal 


Two other eigenvectors are 


(3) (<2) 
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It is often convenient to have small integer values for the components of 
the eigenvectors. To choose an eigenvector with integer components, start 
by writing the equation of the eigenline in the form y = ma. Then choose 
a small non-zero integer value for x such that mz is also an integer. 


To choose such a ‘nice’ eigenvector for the eigenline y = Sy, we need to 
choose a small non-zero integer value for x such that 3x is also an integer. 
In this case, any even integer except 0, will do. 


; 2 
By choosing x = 2, we obtain the eigenvector & ). 


We usually do try to pick small integer values for the components of 
eigenvectors, but in fact any vector of the form (_°). where c € R and 
c # 0, is an eigenvector corresponding to the eigenline y = —a. Similarly, 
any vector of the form ( s.) where c € R and c ¥ 0, is an eigenvector 

2 


corresponding to the eigenline y = 32. 


Activity 2.3 Finding eigenlines 


For each of the following matrices, use the eigenvalues you found in 
Activity 2.2 to find the equations of the corresponding eigenlines. For each 
eigenline, give two eigenvectors. 


@ (Fo) © (02 os) 


Solutions are given on page 53. 


Here is a summary of the process for finding eigenvalues, eigenlines and 
eigenvectors in the form of a strategy. 


Strategy 
To find the eigenvalues, eigenlines and eigenvectors for the 2 x 2 


. a b 
matrix A = (< a 


1. Solve the characteristic equation k? — (a + d)k + ad — bc = 0 to 
find any eigenvalues k. If there are no real solutions to this 
equation, then there are no eigenvalues. 


2. For each eigenvalue k found in Step 1: 
(a) substitute the eigenvalue k into the eigenvector equation 


Ax = kx, where x is the vector (5): 


(b) use the simultaneous equations given by the eigenvector 
equation to find the equation of the eigenline; 

(c) choose a convenient non-zero vector on the eigenline as an 
eigenvector. 
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Remarks 


1. As you will see in Subsection 2.2, it is possible to write down the 
eigenvalues of certain types of matrices without having to solve the 
characteristic equation. 


2. If the eigenline has equation y = mz, then any non-zero vector of the 
form ( s ) is an eigenvector. 
mc 


3. You may be aware that when a quadratic equation has no real 
solutions, it does have solutions that are complex. In some contexts, Complex numbers are studied 
complex eigenvalues are very important. For this course, however, we in Chapter D1. 
are interested only in real eigenvalues; so if there are none, then we say 
that ‘the matrix A has no eigenvalues’ (rather than ‘the matrix A has 
complex eigenvalues’). This corresponds to the definition given on 
page 17, which requires that an eigenvalue be a real number. 


Now try working through the complete process of finding eigenvalues, 
eigenlines and eigenvectors. 


Activity 2.4 Finding eigenvalues, eigenlines and eigenvectors 


Let A be the matrix 


—2 3 
—2 5)° 
Find the eigenvalues and eigenlines of A and, for each eigenvalue, give one 


eigenvector. 


A solution is given on page 53. 


2.2 Special cases of eigenvalues 


In this subsection, we look first at particular types of matrices for which 
the eigenvalues are easy to spot. Then we look at matrices which have one 
eigenvalue equal to zero. Finally, we look at examples of matrices which 
have no eigenvalues or only one eigenvalue. 


Diagonal and triangular matrices 


Activity 2.5 Eigenvalues of diagonal matrices 


(a) Use the characteristic equation to find the eigenvalues of the diagonal 


matrix A = & Y 


Diagonal matrices were 
i: What do you notice? introduced in Chapter B2, 


0 —3 Subsection 2.2. 


(b) Show that the matrix D = ( i has eigenvalues a and d. 


a 
0 d 


Solutions are given on page 54. 
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In Activity 2.5(b), you showed that 2 x 2 diagonal matrices have their 
eigenvalues on the leading diagonal (from top left to bottom right). 


Triangular matrices were In fact, the same argument applies to a triangular matrix. Indeed, a 
int in Chapter B2 : b Sas : 
rr oa ‘ey ee matrix of the form ({ i or (: also has the characteristic equation 


k? —(a+d)k+ad=0, 


which factorises as (k — a)(k — d) = 0. Thus the eigenvalues are again the 
elements a and d on the leading diagonal. 


—2 

0 2 
the eigenvalues are —2 and 2, because this is a diagonal matrix. If we 
recognised this as the matrix of a scaling, then we could argue 
geometrically that the eigenlines would be the x- and y-axes. The next 
activity asks you to verify this algebraically. 


Presented with the matrix ( J, we can now say immediately that 


Activity 2.6 Eigenvectors and eigenlines of a scaling 


Use the eigenvalues —2 and 2 to find the eigenlines of the matrix 


(0 2): 


A solution is given on page 54. 


Zero eigenvalue 


In Chapter B2, Subsection 2.3, you looked at a type of linear 
transformation called a flattening. In particular, you saw that the linear 


. 4 
ut transformation f represented by the matrix A = ( 3) flattens the 


2 3 
plane onto the line 7 — 2y = 0; see Figure 2.4. 
under f 
We now investigate further this flattening by finding the eigenvalues and 


eigenlines of its matrix. 


The matrix A has characteristic equation 


r= Th= 0. 


gv 


ze—2y=0 


This factorises as k(k — 7) = 0, so the eigenvalues of A are 0 and 7. 


Figure 2.4 A flattening onto For eigenvalue 0, the eigenvector equation, 


the line x — 2y = 0 4 6 zr vr 
( 5) =). 


In general, if the determinant . ; ; 

b gives rise to the two equations 4% + 6y = 0 and 2x + 3y = 0. Thus the 
of a matrix 4 eigenline corresponding to the eigenvalue 0 is 2x + 3y = 0. 
representing a linear 


ar For eigenvalue 7, the eigenvector equation, 
transformation is zero, then 


the transformation is a 4 6 t\_7(# 
flattening and the 2 3 y) y )? 
characteristic equation of the : : F 
ennteeas gives rise to the two equations 4% + 6y = 7x and 2% + 3y = 7y. Thus the 
eigenline corresponding to the eigenvalue 7 is x — 2y = 0. 

k? —(a+d)k =0; ; aa ; ; se 
A Since the eigenline 2” + 3y = 0 has eigenvalue 0, every point on this line is 
that is 


scalar multiplied by 0; so every point on this line maps to the origin. 
k(k —(a+d)) =0. 
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You saw in Chapter B2, Subsection 2.3, that this flattening maps every 
point on a line with equation of the form 2x” + 3y = c to the point with 


ois 2 ‘ . : 
position vector ¢ é ). Thus every point of the plane is mapped to a point 


on the eigenline x — 2y = 0. As a check, notice that if @ is the line 

2x + 3y =, which is parallel to the eigenline 2x + 3y = 0, then @ meets 
the eigenline + — 2y = 0 at the point (2c, 4c). This point is indeed mapped 
to the point (2c,c) since the eigenvalue of x — 2y = 0 is 7, as shown in 


Figure 2.5. 
yA yA image 
of € 
f,2n+3y=e¢ 


ay 


image of 
2x + 3y = 0 


Figure 2.5 A flattening onto the eigenline x — 2y = 0, which has eigenvalue 7 


In the above discussion, it was pointed out that every point on the 
eigenline 2x + 3y = 0 maps to the origin. This is an example of an eigenline 
which is not an invariant line, as mentioned in Remark 6 on page 18. 


The next activity concerns another flattening. 


Activity 2.7 Another flattening 


: . 2 
The flattening represented by the matrix A = & ) collapses the plane You saw this in Chapter B2, 


onto the line 3x” — 2y = 0. Activity 2.8. 


(a) Find the eigenvalues of A. 


(b) Deduce that the flattening maps every point of the eigenline 
3x + y = 0 to the origin and that 32 — 2y = 0 is also an eigenline of A. 


Solutions are given on page 54. 


No eigenvalues or one eigenvalue 


The characteristic equation of a 2 x 2 matrix is a quadratic equation, and 
a quadratic equation may have 0, 1 or 2 (real) solutions. You have seen 
several examples of matrices which have two distinct eigenvalues, and 
hence two eigenlines. We now take a closer look at some examples of 
matrices which have either no eigenvalues or only one eigenvalue. Recall 
from Section 1 that most rotations rg have no invariant lines through the 
origin; so we would expect a matrix which represents a rotation (through 
an angle other than a multiple of 7) not to have any eigenvalues. The next 
activity asks you to verify this fact for the rotation through ST. 
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If d= aT, then 


X-( 
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cos 6 
sin 6 


1 


—1 
0 


—sin@ 


). 


cos 0 
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Activity 2.8 A case of no eigenvalues 


ae ‘ . —1 ; 
Use the characteristic equation to show that the matrix Ce 0 ) which 
represents a rotation through an angle of ST, has no eigenvalues. 


A solution is given on page 54. 


We have not yet looked at the case where the characteristic equation has a 
repeated solution. If this happens, then things are a little more 
complicated. When we try to solve the eigenvector equation to locate a 
corresponding eigenline, we may find that there are, in fact, many of them. 


Example 2.2. Many eigenlines 


2 0 
0 2 


(a) Find the eigenvalues and eigenlines of A. 


Consider the matrix A = ( ), which represents a uniform scaling. 


(b) Do your answers agree with your geometric understanding of the linear 
transformation represented by this matrix? 


Solution 


(a) The matrix A is diagonal, so its eigenvalues are the elements on the 
leading diagonal. Thus A has only one eigenvalue, k = 2. The 
eigenvector equation, 


(0 2) (5) =?) 


gives rise to the two equations 27 = 2x7 and 2y = 2y. These equations 
are true for all values of x and y. This implies that every non-zero 
point (x, y) corresponds to an eigenvector, so every line through the 
origin is an eigenline for this matrix! 


(b) We saw in Section 1 that for a uniform scaling, every line through the 
origin is an invariant line, so the above result agrees with our 
geometric understanding of uniform scalings. 


On the other hand, when the characteristic equation has a repeated 
solution there may be only a single eigenline. 


Activity 2.9 A single eigenline 


1 0 


Consider the matrix A = 6 l 


factor 1. 


), which represents a y-shear with 


(a) Find the eigenvalues and eigenlines of A. 


(b) Do your answers agree with your geometric understanding of the linear 
transformation represented by this matrix? 


Solutions are given on page 54. 
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In this course, we consider only 2 x 2 matrices. The definitions of 
eigenvalues and eigenvectors can readily be extended to square matrices of 
larger sizes; but the geometric interpretation of the corresponding 
transformations is harder to visualise. Also, calculation by hand of 
eigenvalues and eigenvectors becomes impractical for large matrices. In 
such cases, a computer can be used. 


Summary of Section 2 


The main theme of this section was how to find eigenvalues and eigenlines 

for 2 x 2 matrices algebraically. In a geometric context, this allowed us to 

identify any invariant lines through the origin for a linear transformation 

represented by a 2 x 2 matrix. This section introduced: 

© eigenvalues, eigenlines and eigenvectors for a 2 x 2 matrix; 

© examples of matrices which have two eigenvalues and hence two 
eigenlines; 

© examples of matrices with no eigenvalues and hence no eigenlines; 

© examples of matrices with one eigenvalue, one example having many 
eigenlines and one example having only one eigenline; 


© the fact that diagonal and triangular matrices have their eigenvalues 
on the leading diagonal. 


Exercises for Section 2 


Exercise 2.1 


Find the eigenvalues (if they exist) and eigenlines of each of the following 
matrices. For each eigenline, give one eigenvector. 


@ (52) © (Co 3356) © ( 4) 
Exercise 2.2 


For each matrix below, identify which type of basic linear transformation it 
represents (rotation, reflection, scaling, shear, flattening) and state how 
many eigenlines it should have, based on your geometric understanding of 
the linear transformation. Verify your answers by finding the eigenvalues 
and eigenlines for each of these matrices. 


ms) (7) (4%) (4) 


all 


3 Using eigenvalues and eigenlines 


Recall that a matrix P is 
invertible if there exists an 
inverse matrix P~! such that 
PP -!=I=P7™!P, where I 
is the identity matrix. 


The matrices in a product 
can be grouped in any 
convenient way provided that 
the order of the matrices is 
not changed. 
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In this section, it is shown that any 2 x 2 matrix A which has two distinct 
eigenvalues can be expressed in the form 


A=PDP", (3.1) 


where D is a diagonal matrix having the eigenvalues of A on its leading 
diagonal, and P is an invertible matrix whose columns are two 
eigenvectors of A. 


As you will see, being able to write A in the form PDP! enables the nth 
power of A to be expressed in the form 

A" =PD"P", 
which provides a systematic method for calculating powers of A. 


Finally we use eigenvalues and eigenvectors to extend our geometric 
understanding of certain linear transformations. 


3.1 Diagonalising a matrix 


In this subsection, you will see how to find a diagonal matrix D and an 
invertible matrix P such that A = PDP™' for any 2 x 2 matrix A which 
has two distinct eigenvalues. 


Consider a general 2 x 2 matrix A and suppose that we can write A as the 
product of the three matrices PDP~', where D is an unknown diagonal 
matrix and P is an unknown invertible 2 x 2 matrix. What do the 
matrices P and D look like? Let 


_ f(a b _(k 0 _ (Pi Pe 
a=(~ i): D=({ i) “o avs 2) 
Since A = PDP“, we can multiply both sides of this equation on the 
right by P to eliminate P~'. We obtain 


AP =PDP 'P 
=PD(P'P) 
=PD, since P"'P =I. 
The equation AP = PD can be written in the form 


(: i) (?: a ees ae a) 
d a 2 qa qe 0 ka)’ 


that is, 
:) & = te ea 
d a 2 kigr  kaq2 ] ° 


( 


We now break up this product by considering each column of P separately. 
The product of A with the first column of P is equal to the first column of 
the matrix PD: 


(oa) i) = (an) 


that is, 


a(?) = ky (7). (3.2) 


io) 


o a 
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The product of A with the second column of P is equal to the second 
column of the matrix PD: 


Ce aa 


that is, 


A (”) = ky ) (3.3) 


Equations (3.2) and (3.3) are exactly of the form of the eigenvector 

equation of A! So these equations are satisfied provided that k, and kz are 

eigenvalues of A, with corresponding eigenvectors G 2 and & : ). 
1 2 

Thus, to obtain the equation AP = PD, we can take D to have the 

eigenvalues of A as its diagonal elements and the columns of the matrix P 

to be corresponding eigenvectors of A. 


To write A in the form PDP™', we need to calculate P~! as well. Since the columns of P are 
vectors in different directions, 
det P is non-zero, so P~+ 


exists. 

Example 3.1 Finding the matrices P and D 

Express the matrix 
1 2 
a=(3 9) 

in the form PDP™', where D is a diagonal matrix. 
Solution 
From the preceding discussion we know that the eigenvalues of A are the 
diagonal elements of D and corresponding eigenvectors of A are the 
columns of P. We found the eigenvalues and corresponding eigenvectors of 
A in Section 2. See pages 18 and 21. 


: ‘ : 2 
© The eigenvalue 4 has corresponding eigenvector & ). 


; : : 1 
© The eigenvalue —1 has corresponding eigenvector & }. 


Thus we can choose 
4 0 2 1 
De, a) and eS ae 


Also, det P = —5, so P~! exists and is given by 


ae go aes 
Ste 
. =3 (53 _ 


_lfl 14 
~~ 5\3 —2)° 


Thus we have A = PDP™. 
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For example, we could have 


4 —2 
wed P= ( > fr 


1 2 2 
. —1 _— 
which P~* = 0 € a) 


The columns of P are 
eigenvectors in different 
directions, so P is invertible. 
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Comment 


You can check the above result by multiplying the matrices 


 f2 T\f4 Oy)171 4 
EDF =F a DEG a 


When there is a number (4 here) involved in a matrix product, it is usually 
simpler to take the number to the front and perform the scalar 
multiplication last. Also, remember that you can multiply a matrix 
product of the form ABC as A(BC) or as (AB)C. Thus, for example, 


(3 -1) (0 -1)5(s 2) =3(2 4) (3 —) 
=F o) 


oft. 2 eee 
=|3 9]}> as required. 


Note that the order in which the eigenvalues are entered in the leading 
diagonal of D matters. If we had written down the matrix D as 


—1 0 
0 4)’ 
then a matrix P such that A = PDP™! would be & - In other 


words, the columns of the matrix P have to be interchanged! The order of 
eigenvalues and eigenvectors must be matched up. 


Note that there are many invertible matrices P which could be used in 
equation (3.1), since there are many possible eigenvectors for each 
eigenvalue. 


We have now developed the following strategy to express a 2 x 2 matrix A 
with two distinct eigenvalues in the form A = PDP™', where D is a 
diagonal matrix. We call this process diagonalising the matrix A. 


Strategy 


To express a 2 x 2 matrix A with two distinct eigenvalues in the form 
PDP, where D is a diagonal matrix. 


1. Find the eigenvalues k, and kz of A. 


2. Define the diagonal matrix D by D = & : ) 

2 
3. Find the eigenlines corresponding to the eigenvalues ky; and kz. 
4. Choose 


© an eigenvector e ; for the eigenvalue k; 
1 


© an eigenvector @ . for the eigenvalue ky. 
2 
5. Use these eigenvectors to define the matrix P by P = @ i ). 
1 


6. Calculate the matrix P~. 
Then A= PDP". 


SECTION 3 USING EIGENVALUES AND EIGENLINES 


Remarks 


1. When applying this strategy, remember that you must match the order 
of the eigenvalues on the diagonal of D and that of the corresponding 
eigenvectors in the columns of P. It does not matter which order you 
take, as long as you are consistent. 

2. As a convention in this chapter, the eigenvalue of larger magnitude will 
be written in the first column of D; so if the matrix A has eigenvalues 


k, and ky with |ki| > |k2|, then D = & k ) 
2 


3. There are many choices of eigenvectors for the columns of the matrix 
P. In order to keep the arithmetic simple, try to choose eigenvectors 
that have small integer components. 


4. It is important to remember that this strategy works only for matrices 
A which have two distinct eigenvalues, so you cannot apply it to 
matrices which represent linear transformations such as rotations or 
shears. 


Activity 3.1 Using the strategy 


Express the matrix A = G : 


) in the form PDP™', where D is a 
diagonal matrix. 

A solution is given on page 55. 

Comment 


Note that in the solution, we could have swapped the order of the 
eigenvalues and used the matrices D and P, where 


—2 0 2 3 
p=( 0 ) and ee ae 


We could also have chosen different eigenvectors for the columns of P. 


3.2 Matrix powers 


Now that we know how to diagonalise a 2 x 2 matrix having two distinct 
eigenvalues, we look at why diagonalisation is of use in calculating matrix 
powers. The following activity illustrates a useful property of diagonal 
matrices. 


Activity 3.2 Powers of a diagonal matrix 


Let A= G ) and D = & Eat Calculate the matrices 
A’, A?,D? and D’. 


Solutions are given on page 55. 


el 


Remember not to change the 
order of the matrices when 
regrouping. 
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As you will have noticed, calculating powers of a general 2 x 2 matrix can 
be tiresome, whereas calculating powers of a diagonal matrix is easy. In 


: a we have 


general, for a diagonal matrix D = é d 


n_ [a 0 
ie ( 0 - ; 
Now suppose that we have a 2 x 2 matrix A which can be written in the 
form A = PDP“, where D is a diagonal matrix. Then 


SAA 
=(PDP™')(PDP™'), since A= PDP", 
=(PD)(P“'P)(DP™'), regrouping the matrices, 
=(PD)I(DP™"), since P-'P =I, 
=PD'P~. 

Similarly, 

A’ = AA? 
=(PDP')(PD’P™'), from above result for A’, 
= (PD)(P7'P)(D’P“') 
= (PD)I(D’P“') 
=PD'P*", 


In general, we have the following result. 


Calculating powers 
If a 2 x 2 matrix A can be written in the form 


A=PDP™', where D is a diagonal 2 x 2 matrix, 
then 
A" =PD"P™, forn=1,2,3,.... 


Writing a matrix A in the form PDP! where D is a diagonal matrix 
greatly simplifies the calculation of powers of A. To calculate A”, we first 
calculate D”, which is straightforward since D is a diagonal matrix, and 
then multiply the three matrices P, D” and P™!. 


Example 3.2 Matrix powers 


As we saw in Example 3.1, the matrix A = € _ can be written in the 


form PDP! where 


_f2 4 (4 0 » ft 4 
P=(5 4) B=, i) and P =3(5 2) 


Use this information to calculate A®, without multiplying powers of A 
directly. 
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Solution 


First we calculate D®: 


p= (5 (a5) 
al ae 


Then we calculate A® using the fact that A” = PD"P7!: 
=P Pp 

[2 1 4096 0 i, 1 1 

~ \3 -1 0 1/5\3 -2 

a 8192 1 1 1 

~ 512288 -1/\3 —-2 

= I 8195 8190 

— 12285 12290 


5 


1639 1638 
2457 2458 } ° 


Now it is your turn to calculate a matrix power. 


Activity 3.3 Matrix powers 


Consider the matrices 
Re = jee € a) and D = € oo. 
(a) Calculate the matrix P™. 
(b) Verify that 
A=PDP™,. 
(c) Hence calculate A*, without multiplying powers of A directly. 


Solutions are given on page 55. 


When calculating a matrix power in this way, it may seem that having first 
to calculate the matrices P, D and P™! is a lengthy preliminary 
procedure. There are two reasons why diagonalisation is a useful way to 
calculate matrix powers by hand: speed and accuracy. 


If n is a large number, then calculating D” is much faster than calculating 
A”. Admittedly, diagonalisation also involves having to first calculate the 
matrices P, P~! and D, as well as multiplying three matrices together at 
the end, but if n is a large number, then diagonalisation is still a much 
quicker way of calculating matrix powers than direct computation. 


oS 
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For example, with A, P and D as in Example 3.2, 


pw (4° 0) _ (1048576 0 
~Lo (pe) > - i 


Then, using exactly the same form of calculation as for A® in Example 3.2, 
we obtain 
419431 419 430 
10 _ 10p-1 _ 
Skee (oe Sv ; 


Calculating A?!° directly would involve many more multiplications and 
additions, so it offers more potential for making arithmetical errors. 


Another point to consider here is accuracy. Every time a calculation is 
carried out by a calculator as an adjunct to hand calculation, the accuracy 
of the answer is limited by the number of digits stored in memory. The 
more calculations are done, the more errors creep in due to rounding. 
Calculation of matrix powers using diagonalisation involves fewer 
calculations and hence gives more accurate answers. 


3.3 Using eigenvalues and eigenlines geometrically 


In Chapter B2, you have seen how a linear transformation may be 
interpreted geometrically by sketching the image of the unit grid. Another 
grid, based on eigenlines, provides an alternative way to interpret a linear 
transformation. 


Here we consider a linear transformation f represented by a 2 x 2 

matrix A that has two distinct non-zero eigenvalues. We already know 
that the image of a point on an eigenline of A lies on that eigenline, being 
scaled by the corresponding eigenvalue. But what happens to points that 
are not on an eigenline? 


For example, consider the linear transformation f represented by the 


matrix A = i a which, as we have seen in Activity 3.1, has 
eigenvalues k, = 3 and ky = —2, and corresponding eigenlines @, and 4 
with equations y = 5x and y = —32, respectively. 


On the left of Figure 3.1 the eigenlines of A are shown, in the domain of f. 
Let P be a point in the domain of f that does not lie on an eigenline. The 
following construction shows what happens to P under f. 


Draw a parallelogram, based at the origin, by joining P to each eigenline 
by a line segment parallel to the other eigenline, as shown. Let the vertices 
of this parallelogram be O, A (on £2), B (on ¢,) and P, as shown. Each of 
the sides OA and OB of this parallelogram OAPB is scaled under f along 
the eigenline on which it lies by the corresponding eigenvalue. Since linear 
transformations preserve parallelism, the image of OAPB is a suitably 
scaled parallelogram. Let A’ and B’ denote the images of A and B, 
respectively. To construct this image parallelogram, mark A’ on the line ¢5 
twice as far from O as A but on the opposite side of line 4; to A, and mark 
B’ on the line ¢; three times as far from O as B and on the same side of 
line £ as B, as shown on the right of the figure. (In terms of eigenvectors, 


—_ —_> = = 
OA' = —20A and OB' = 30B.) Now draw the parallelogram with sides 
OA’ and OB’. The fourth vertex of this parallelogram P’ is the image of P. 
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Figure 3.1 Constructing P’ 


We have located P’ by making use of the eigenvalues and eigenlines of A. 
The effect of f on P can be described as follows: P is scaled by the 
factor 3 in the direction of the eigenline @; and by the factor —2 in the 
direction of the eigenline 2. In this sense, the linear transformation f can 
be described as a generalised scaling. A linear transformation represented 
by a 2 x 2 matrix that has two distinct non-zero eigenvalues (and hence 
two distinct eigenlines) is called a generalised scaling. 


The above construction applies to any point P. In particular, consider P If P lies on the eigenline £1, 
to be such that OA and OB are of unit length. Then the parallelogram say, then the construction 

OAPB can be thought of as a ‘unit parallelogram’, which can be used as__ yields the line segment OP’, 
the basis of a ‘unit grid’ of parallel lines, each of which is parallel to an along £;. The scaling in the 


direction of £2 has no effect 


eigenline, as shown on the left in Figure 3.2. The image of this grid under ‘ 
on P. 


f is shown on the right: each point of the original unit grid and each 
parallelogram has undergone the same generalised scaling. 


ya yA 


Rv 


Figure 3.2 A grid of lines parallel to the eigenlines of A, and its image 


In the next activity, you are asked to apply the above construction. 


oo 


Zero eigenvalues are excluded 
because they correspond to 
flattenings. 


If P lies on £9, then P’ is also 
on fy. So there is a sense in 
which P’ lies on the same side 
of é) as P. 
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Activity 3.4 Constructing and confirming 


—2 
—2 5 
eigenlines y = 2a and y = a as shown in Figure 3.3. 


(a) Construct on Figure 3.3 the image P’ of the point P(1,1) under the 
linear transformation represented by A. 


The matrix A = ( ) has eigenvalues 4 and —1, with corresponding 


(b) Check your construction by calculating P’. 


Figure 3.3 Locating P’ 


Solutions are given on page 55. 


This subsection concludes with some properties of generalised scalings, and 
an activity in which you will use the properties. 


Properties of generalised scalings 
Let the linear transformation f be represented by a 2 x 2 matrix A 
that has two distinct non-zero eigenvalues k, and kz (that is, f is a 
generalised scaling) with corresponding eigenlines ¢; and £2. Let P be 
a point of R® with image P’ under f. 
(a) (i) If k, > 0, then P’ lies on the same side of f as P. 

(ii) If k; < 0, then P’ lies on the opposite side of £2 to P. 


(b) The distance from P’ to £y is |k,| times the distance from P to é. 


Remarks 


1. Corresponding properties are obtained if k; and @, are replaced by kz 
and ¢. 

2. In property (b), the distance from a point to the line @, may be taken 
to be the perpendicular distance from the point to the line, or to be 
the distance in the direction of the line ¢, from the point to the line £5. 
See Figure 3.4. 
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(a) perpendicular distance (b) distance in direction of & 
Figure 3.4 Distance from a point to 02 


3. If|ki| >1, then P’ is further from @, than P is. See Figure 3.5(a). 
If |k,| < 1, then P is further from @, than P’ is. See Figure 3.5(b). 


ya bo 


(a) [Ai] > 1 (b) |i] <1 


Figure 3.5 Comparing (perpendicular) distances from £3 


Activity 3.5 Locating images 


This activity concerns the same matrix and linear transformation as in 
Activity 3.4. 


(a) Use properties of generalised scalings to describe the location of 
(i) the image P’ of the point P shown in Figure 3.3; 
(ii) the image P” of the point P’. 

(b) Calculate P”. 


Solutions are given on page 56. 


ral 
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Summary of Section 3 


This section has introduced: 


© the process of diagonalisation — if a 2 x 2 matrix A has two distinct 
eigenvalues, then it can be written in the form A = PDP’, where D 
is a diagonal matrix whose elements are the eigenvalues of A and P is 
an invertible matrix whose columns are corresponding eigenvectors 


of A; 


© the use of diagonalisation to compute powers of A quickly and 
accurately, by means of the result 


A®=PD°P™, fora#=1,3.3,.3< 


© the notion of a generalised scaling and its properties. 


Exercises for Section 3 


Exercise 3.1 


(a) Express each of the following matrices in the form PDP™', where D is 
a diagonal matrix. (You can use the solution to Exercise 2.1 here.) 


0 (3) 2S 290) aw (2 2) 


(b) For the matrix in part (a)(i), check your calculation by multiplying 
out the matrix product PDP™' that you have found. 

(c) For the matrix A in part (a)(i), find A° without multiplying powers of 
A directly. 


Exercise 3.2 


: : 1 ; . : 
Consider the matrix A = & ) whose eigenvalues and eigenlines you 


5 2 
calculated in Exercise 2.1(a). 


(a) Use properties of generalised scalings to describe the location of the 
image P’ of the point P(1,—2) under the linear transformation 
represented by A. 

(b) Calculate P’. Draw a sketch to show that the locations of P and P’ 
match your description in part (a). 


4 Iterating linear transformations 


So far in this chapter, we have looked at the effect of a single application of 
a linear transformation f. But what happens to the points of R’ if we 
repeatedly apply f; that is, if we apply the process of iteration? 


In Chapter B1, you saw iteration sequences in the context of real 
functions. Here we look at iteration sequences for linear transformations. 
Given a linear transformation represented by a matrix A and an initial 
point (2, yo), represented by the vector xo, the iteration sequence 
generated by A with initial point (29, yo) is the sequence x,, of points 
represented by the vectors 


Xo, X,=AXp, X2=AX), X3 = AX, 


This is the sequence of points in R° obtained by repeated application of 
the linear transformation represented by the matrix A, starting with the 
initial point (29, yo). For example, if the matrix A represents the 

1 Baa 
0): then the sequence x, appears as in Figure 4.1. 


rotation rg and xg = 


Figure 4.1 The sequence x, 
Given an initial point represented by the vector xo, the corresponding 
iteration sequence generated by the matrix A has recurrence relation 
CoS ae, Cah 
This sequence can be expressed in closed form as 
= Ake (= 1,2,3)<55): 


In this section we consider the behaviour of iteration sequences generated 
by matrices which have two distinct eigenvalues. We first look at such 
iteration sequences in which the initial point is on one of the eigenlines 
of A. 


x2 => A?’xo, 


x3 = A®xg,.... 


Bg 


In particular, see page 21. 


f is a generalised scaling. 
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Activity 4.1 Iteration of a point on an eigenline 


Consider the matrix 


a=(53) 


Earlier, you saw that A has eigenvalues 4 and —1 with corresponding 
eigenlines y = 2a and y = —a. Let xo be the vector representing the point 
(2,3), which lies on the eigenline y = 3a. Calculate the first three points of 
the iteration sequence generated by A with initial point (2,3). 


A solution is given on page 56. 


Now consider a general 2 x 2 matrix A having two distinct eigenvalues, k, 
and kz, with corresponding eigenlines ¢; and @5, respectively. If we start 
with a point (Zo, yo) on the eigenline @, and apply the linear transformation 
f represented by A, then the resulting point is (k,29, iyo), which is also 
on the eigenline ¢,. If we apply the same linear transformation to the point 
(k129, kyyo), then we obtain the point (k729,k?yo). If we continue iterating 
this linear transformation, then after n applications of the transformation 
we obtain the point (k?29, k7'yo). Similarly, if we start with a point (2, y) 
on the eigenline 2, then after n iterations, we reach the point (k}ap, ky yp). 
These sequences are illustrated in Figure 4.2, for the case k, > 1, ky < —1. 


YA 


(k326, kup) 
(Kio, Kiyo) 


(Kizo, Kiyo) 
(Kio, kyo) 
(wo, yo) 


Ry 


(0, vo) 


(k326, k3y6) 


eigenline (2 
with 
eigenvalue 
kg < -1 


eigenline @1 
with 
eigenvalue 
ky > 1 


Figure 4.2 Iteration sequences on eigenlines 
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The table below describes the long-term behaviour of iteration sequences 
x, generated by the matrix A of a generalised scaling. Here the non-zero 
initial point x, lies on an eigenline ¢ corresponding to a (non-zero) 
eigenvalue k of A. The table is based on the properties of generalised 
scalings given in Subsection 3.3. 


k Long-term behaviour of x, on 
k>1 Xn moves away from (0,0), on the same half of & as xo 
k= Xn = Xo, for n = 0,1,2,... (a constant sequence) 


A line through the origin is 
divided into two halves by the 
origin. 


0<k<1 x, moves towards (0,0), on the same half of @ as xo 
-1<k<0O x, moves towards (0,0), alternating between the halves of ¢ 
k=-1 Xp, alternates between + xo 

k<-l X, moves away from (0,0), alternating between the halves of ¢ 


In the next activity you will use this table. 


Activity 4.2 Long-term behaviour 


For the matrix A of Activity 4.1, use the table above to describe the 
long-term behaviour of the iteration sequences generated by A with the 
following initial points. 

(a) (2,3) — (b) (2, -2) 


Solutions are given on page 56. 


If we start with an initial point on an eigenline, then the corresponding 
iteration sequence consists only of points which are on that eigenline. But 
what can we say about an iteration sequence whose initial point is not on 
an eigenline? Here is an example. 


Example 4.1 Iteration of a point not on an eigenline 


Consider again the matrix A = é >) and the point (2,1) represented Note that the point (2,1) 
satisfies neither of the 
by the vector xp = (; ). eigenline equations y = —x 
1 and y = 3a. 


(a) Calculate the first four points in the iteration sequence 
Gag SAR, HH 0,1,2...;5) 
with initial point (2,1). 


(b) Find a formula in terms of n for the vector x, which represents the 
(n+ 1)th point in this iteration sequence. Use that formula to 
calculate the 10th point, represented by the vector xg. 
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yA yA 
(76, 116) 
40- 
© (20, 28) 
207 
(4,8) 
2,1 
e251) | | 
20 49 © 


Figure 4.3 The first four 
points 


We know these matrices from 
Example 3.2, in which A® 
was calculated. 


In these matrices, to save 
space, products like 2 x 4” 
have been written as 2(4)”. 


Remember that the (n + 1)th 
point of this sequence is 
represented by the vector xp. 
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Solution 


(a) The first point in this iteration sequence is (2,1). 


The second point is represented by the vector 


n-tn-(3 2)(2)-(9) 


The third point is represented by the vector 


nn an=(2 3)()- (2) 

The fourth point is represented by the vector 
ty & 3) & ~ (is) 

Thus the first four points in this iteration sequence are 
(2,1), (4,8), (20, 28), (76, 116); 

see Figure 4.3. 


(b) To find a formula in terms of n for x,,, we use the closed form 


x,, = A”x, and the fact that A = PDP™', where 


_f2? 1 (4 0 leit 4 
Pag 2) DS(5 7) MPS, Sl: 


We first calculate a formula for A”: 


A” =PD"P™ 


=(5 a) (> cay) (3) 


(-1 
1 /2(4)" —1)" 1 1 
(sue cr) (¢ 2) 
1 Ga +3(-1)” 2(4)"— ae 
5 \3(4)" — 3(-1)"_ 3(4)" + 2(-1)” 
Now we calculate x, = A”xo: 
. 1 /2(4)"+3(-1)” 2(4)” — 2(-1)" 2 
= Arm = 5 (Saye ate ate ty») (2) 
_ te + 6(—1)" + 2(4)" — a) 
5 \ 6(4)" — 6(—1)" + 3(4)" + 2(-1)” 
1 Gay ee 
5 \9(4)” — 4(-1)” 


To find the 10th point in this iteration sequence, we substitute n = 9 
into the expression for x,: 

1 6(4)° + 4(-1)9 
5 \ 9(4)° — 4(-1)° 
1 ( 1572 a) 


~ 5 \ 2359300 


_ (314572 
~ \ 471860 ) ° 


Xo = 
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Comment 

A quicker way to calculate x,, = A”xo is as follows. 
A”x) = PD"’P!x, 

= (PD")(P“'xp 


ae, 
of 
a 

| 

HK © 
ea 

3 
Se 
—l 
ES | 
a 
wre 
| 

| Oe 
Se 
oS 
re bo 
4 
_— | 


( 
= 5 (ou) aca): 


A disadvantage is that PD"”P~ is not calculated explicitly, and so is not 
available for future use; see Activity 4.3. 


From Figure 4.3, it appears that the points of this iteration sequence are 
moving up and to the right, away from the origin, in roughly a straight 
line. The ratio y,,/2, for each of the calculated points (x, y,) of this 
iteration sequence gives us an indication of the direction in which the 
iteration sequence is moving. 


Point (an, Ym) | (2,1) (4,8) (20,28) (76,116)... (314572,471860) 
Ratio yn/an | 0.5 2 14 1.526316 ... 1.500 006 


It appears that the ratio y,/x» is tending towards the value 1.5, so we 
might conjecture that in the long term, the iteration sequence is moving 
up and to the right, away from the origin, in the direction of the straight 
line through the origin with gradient 1.5, that is, the eigenline y = 20. 


Is this behaviour repeated if the initial point is changed to another point 
1 2 


3 9 ) ? In the next activity, you 


not on either of the eigenlines of A = ( 


are asked to investigate this question. 


Activity 4.3 A different initial point 


Consider the matrix A from Example 4.1, and the initial point (—3, —1). 


(a) Calculate the first four points in the iteration sequence 


xii Ax, (= 0,1,2,4.5) 


with xp = ei 


(b) Use the matrix A” calculated in Example 4.1 to find a formula in 
terms of n for the vector x/,, which represents the (n + 1)th point in 
this iteration sequence. Use that formula to calculate the 12th point in 
this iteration sequence. 


(c) Calculate the ratio y/ /x', for each of the points (z/,, y/,) found in 
parts (a) and (b). Do these ratios appear to tend to a particular value? 


Solutions are given on page 56. 
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The first four points for each of the two iteration sequences x, and xi, 
calculated above, are plotted in Figure 4.4. 


Rv 


i 
x3 40 
(=101,,—155) 


Figure 4.4 Two iteration sequences 


In Figure 4.4, it appears that both iteration sequences are tending in the 
direction of the eigenline y = oo. (Note that ‘tending in the direction of 
the eigenline y = 3 3 does not imply that the sequence is ‘tending to the 
eigenline’. By property (b) of generalised scalings, the eigenvalue —1 
associated with the eigenline y = —x implies that each term of the 
sequence is the same distance from the eigenline y = 3a, and consecutive 
terms lie on opposite sides of that line.) 


We now show this conjecture to be true for the iteration sequence x,. In 
Example 4.1, we found a formula for the vector x,,, namely 


Cin) =3 (ote ap): 


Using this formula, we obtain the ratio 
Yn  9(4)" —4(-1)” 


fm, 6(4)" +4(—1)”" 
The value of —4(—1)” is much smaller than 9(4)" and the value of 4(—1)” 
is much smaller than 6(4)", so in a sense the terms involving 4” ‘dominate’ 
for large values of n. To be precise, we divide both numerator and 
denominator by 4”, to obtain 


Ln 6 + 4(—4) 


n 
nm 
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As n gets large, (—4)” tends to zero, so The sequence (—+)” is a 
reometric sequence with 

Yr 9-4(0) 3 . a 

— as n— Co. common ratio q: 


tn,  6+4(0) 2 


Since the ratio y,,/x, tends to 1.5 (the gradient of the eigenline y = 32), 
the iteration sequence with initial point (2,1) tends in the direction of this 
eigenline, as conjectured. 


In fact, any iteration sequence generated by the matrix 


a=(9 3) 


which starts at a point not on an eigenline tends in the direction of the 
eigenline y = Sy. However, the algebra required to prove this result is 
quite tedious, so instead a geometric explanation of this long-term 
behaviour is given. 


Recall, from Subsection 3.3, that a linear transformation represented by a 

matrix which has two distinct non-zero eigenvalues is a generalised scaling. 
In this particular case, each time the linear transformation represented by 

the matrix A is applied, the points of the plane are scaled by the factor 4 

in the direction of the line y = 3a and scaled by the factor —1 in the 


direction of the line y = —x. So successive points of the iteration sequence 
are scaled by the factor 4 in the direction of the line y = 32, and by the 
factor —1 in the direction of the line y = —2. 


Using property (a) of generalised scalings, we deduce that successive points See page 36. 
of an iteration sequence 


© stay on the same side of the line y = —z, since 4 > 0; 
© move from one side of the line y = $x to the other, since —1 < 0. 


Figure 4.4 shows the two iteration sequences x,, and x’, exhibiting this 
behaviour. 


The following result enables the long-term behaviour of an iteration 
sequence generated by a matrix representing a generalised scaling to be 
described. Part (a) is a restatement of property (a) of generalised scalings 
for the terms of an iteration sequence. Part (b) contains criteria for 
deciding whether an iteration sequence moves away from (0,0), as in the 
example discussed above, or whether it moves towards (0,0). Part (c) 
describes the long-term behaviour of the ratio y,,/x,, in certain cases (the 
behaviour of this ratio was calculated in the example above). 


In this result, the notation max{a, b} is used to denote the maximum of 
two real numbers a and 6b. For example, 


max{2,-5}=2 and max{l,1}=1. 
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Iteration properties of generalised scalings 


Let the linear transformation f be represented by a 2 x 2 matrix A 
that has two distinct non-zero eigenvalues k, and ky with 
corresponding eigenlines @; and @. Let (x9, yo) be a point of R? which 
is not on an eigenline of A, and let (x, y,) be an iteration sequence 
generated by A, with initial point (29, yo). 
(a) (i) If k, > 0, then all the points of (x,,,y,) lie on the same side of 
fy as (Zo, Yo). 
(ii) If k, < 0, then the points of (7, y,) alternate between 
opposite sides of 05. 


(b) (i) If max{|k,], |k2|} > 1, then the sequence moves away from 


? 


(ii) If max{|k1|, |ko|} <1, then the sequence moves towards (0,0). 
(c) If |kiy| > lka|, then 


If the eigenline ¢; is x = 0, ue > masn— oo, 

then Tn 
Yn where m is the gradient of ¢;. 
— +o as n—- oo 
Xn 


Remarks 

1. Properties corresponding to (a) and (c) are obtained if k, and é are 
replaced by ky and ¢,, and vice versa. 

2. If |ki| > |k2|, then k, is called the dominant eigenvalue and /; is 
called the dominant eigenline. Property (c) is called the Dominant 
Eigenvalue Property. 


3. The long-term behaviour specified in property (c) is often expressed as 
‘the sequence (2p, Yn) tends in the direction of the line ¢,’. 


4. Note that the case max{|k,|,|k2|} = 1 is not covered by property (b), 
and the case |k,| = |k2| is not covered by property (c). 


In this section you will consider only matrices for which no eigenvalue has 
modulus less than 1. The next activity concerns such a case. (In Section 5, 
you will apply the properties to other cases.) 


Activity 4.4 An iteration sequence 


Consider the linear transformation represented by the matrix 


—2 1 
dee: ( : 1) 
This matrix has eigenvalues —3 and 2, with corresponding eigenlines 


y = —x and y = 42. (You are not asked to establish these facts.) 


Let (2,1) be the initial point of an iteration sequence (2,,, y,) generated by 
A. Note that (2,1) does not lie on either of the eigenlines of A. 


Use iteration properties (a), (b) and (c) of generalised scalings to describe 
the long-term behaviour of the sequence (2, Yn). 


A solution is given on page 57. 
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An application to population models 


This section concludes with a brief description of an application of matrix 
iteration to the long-term structure of subpopulations in population 
models. For example, a particular model of the UK population involves 
just two subpopulations — juveniles (those less than 15) and adults (the 
rest). Let J, denote the number of juveniles and A, denote the number of 
adults at the beginning of year n. 


The relationship between these subpopulations in successive years is 
modelled by the recurrence relation 


Insi \ — (0.9326J, + 0.01724, : 
Cr) ~ Ge + 0.98644, (n = 0,1,2,...). 


This equation can be expressed in matrix form as 


Phi = Mp,,, 


J _ {0.9326 0.0172 — awe 
where p,, = ce and M = Cie ao Thus p,, = M"p,, where 


Pp represents (Jo, Ao). 
It turns out that the matrix M has distinct eigenvalues: 


k, = 1.0027, witheigenline y= 4.07752; 
ky = 0.91627, with eigenline y= —0.949 62x. 


(The values given are correct to five significant figures.) 


Since k, > ky > 0, the dominant eigenvalue is k,. According to the 
Dominant Eigenvalue Property, we have (for most initial subpopulations) 


—" _, 4.0775 as n > oo. 


So, according to the model, in the long-term there will be approximately 
four times as many adults as juveniles in the UK population. 


This example is just one illustration of the many applications of matrices 
to the study of population models. 


Now would be a good time to watch the optional band on DVDO00095. 
Now watch band B(iv), ‘Weaving spirals’. 


The Video Band shows Mathcad being used. This is a previous version of 
Mathcad from the one currently used in MS221, though there is no 
difference in the way the two versions work, as far as this Video Band is 
concerned. 


The rest of this section will 
not be assessed. 


This population model, with 
n being the number of years 
after June 1990, was 
considered in MST121 
Chapter B2, Section 3. 


This limit agrees with the 
result of numerical calculation 
given in MST121 Chapter B2. 
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Summary of Section 4 


In this section we investigated iteration sequences for linear 
transformations represented by 2 x 2 matrices having two distinct 
eigenvalues. 


This section has introduced: 
© iteration sequences generated by a 2 x 2 matrix A 
Sep Ae, (=H 01, ei ccc) 
with initial point (xo, yo) represented by xo; 


© the use of diagonalisation of A (A = PDP“, where A has two 
distinct eigenvalues) to express the closed form x, = A”X of an 
iteration sequence in terms of n; 


© results describing the long-term behaviour of iteration sequences 
generated by matrices representing generalised scalings. 


Exercises for Section 4 


Exercise 4.1 


Determine the matrix A”, where 


a=(i 6) 


has diagonalisation 


a=(1 a) (0 2) 


Exercise 4.2 


(i 4). 


) whose eigenvalues and eigenlines you 


1 
5 


6 
5 2 
calculated in Exercise 2.1(a). 


(a) Which of the points (1,3), (2,2) and (1, —5) lie on an eigenline of A? 


(b) Describe the long-term behaviour of iteration sequences (2n, Yn) 
generated by A with each of the following initial points. 


(i) (1,3) 
(ii) (2,2) 
(iii) (1, —5) 


Consider the matrix A = ( 


5 Iterating linear transformations with the 
computer 


In this section you will need computer access and Computer Book B. 
oO 


In Section 4, we looked at the long-term behaviour of iteration sequences 
for linear transformations which have two distinct non-zero eigenvalues. In 
this section, we consider such long-term behaviour for more general linear 
transformations. We also look at the iteration of some affine 
transformations. Recall from Chapter B2 that an affine transformation is a 
transformation of the form 


f:R —R 
xh Ax+a, 
where A is a 2 X 2 matrix and a is a vector. 


We use the computer to investigate the iteration of linear and affine 
transformations, and discover how iteration can lead to pictures of natural 
objects such as ferns. 


Refer to Computer Book B for the work in this section. 


Summary of Section 5 


This section has used the computer to help discover patterns occurring in 
the iteration of linear and affine transformations. 
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In this chapter, you met eigenvalues, eigenlines and eigenvectors of a 2 x 2 
matrix. 


You saw that a matrix A which has two distinct eigenvalues can be 
diagonalised. The diagonal form enables powers of A to be calculated 
easily. 


The long-term behaviour of iteration sequences generated by certain linear 
transformations was described. 


Learning outcomes 


You have been working towards the following learning outcomes. 


Terms to know and use 


Fixed point, invariant line, eigenvalue, eigenvector, eigenline, 
characteristic equation, eigenvector equation, diagonalising a matrix, 
generalised scaling, dominant eigenvalue, dominant eigenline. 


Notation to know and use 


max{a, b} 


Mathematical skills 


© Determine the fixed points and invariant lines (through the origin) of 
certain linear transformations. 


© Determine the eigenvalues of a 2 x 2 matrix using the characteristic 
equation. 


© Find the eigenline corresponding to a given eigenvalue, and choose 
‘nice’ eigenvectors. 


© Write down the eigenvalues of diagonal and triangular matrices. 


© Use the strategy for diagonalising a 2 x 2 matrix A which has two 
distinct eigenvalues, and hence write the matrix A in the form 
PDP, where D is a diagonal matrix. 


© Compute powers of a 2 x 2 matrix which is expressible in the form 
PDP, where D is a diagonal matrix. 


© Given a 2 x 2 matrix A with two distinct eigenvalues and an initial 
point (29, yo), describe the long-term behaviour of the iteration 
sequence generated by A. 


Mathematical awareness 


© Know that 2 x 2 matrices might have no eigenlines, one eigenline, two 
eigenlines or many eigenlines. 


Mathcad skills 


© Investigate iteration sequences of linear and affine transformations. 


SUMMARY OF CHAPTER B3 


Summary of Block B 


This block has introduced several key topics: 


© iteration sequences, with recurrence relations of the form 


C4 = {ie C=O ees) 
where f is a real function and 29 is a given initial term; 
© the Binomial Theorem; 


© linear transformations f : R? —> R° whose rules are of the form 
f(x) = Ax, where A is a 2 x 2 matrix; 


© affine transformations f : R? —> R° whose rules are of the form 
f(x) = Ax +a, where A is a 2 x 2 matrix and a is a vector; 


© fixed points and invariant lines for linear transformations; 


© 


eigenvalues, eigenlines and eigenvectors for 2 x 2 matrices; 


. . . . 2 . ‘ 
© iteration sequences of points in R”, with recurrence relations of the 
form 


X09 =F (eq) = 0, 122,222), 

where f is a linear transformation or an affine transformation, and xg 

is a given initial point. 
In studying the long-term behaviour of iteration sequences generated by 
real functions, you saw the importance of the fixed points and p-cycles of 
the function f being iterated, and how these are classified using the 
gradient of the graph of f. The gradient of the graph of a real function is 
the fundamental concept in calculus, and you will learn much about this 
topic in Block C. 


In studying linear transformations, you saw that the geometric nature of a 
linear transformation f can be related to properties of the matrix A that 
represents f. This relationship enabled us, for example, to describe the 
long-term behaviour of iteration sequences of points generated by certain 
linear transformations known as generalised scalings. You will see more 
about this relationship between linear transformations and their matrices 
in Block D. 
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Solutions to Activities 


Solution 1.1 


(a) The rotation r3,/2 has only one fixed point, the 
origin. 

(b) The reflection gz has a line of fixed points, the 
y-axis. 

(c) The scaling with factors 1 and —1 has a line of 
fixed points, the z-axis. 


(d) The y-shear with factor 4 has a line of fixed 
points, the y-axis. 


Solution 1.2 


(a) The rotation r3,/2 has no invariant lines through 
the origin. 


(b) The reflection q,/2 has two invariant lines 
through the origin, the y-axis (which is a line of 
fixed points) and the z-axis. 

(c) The scaling with factors 1 and —1 has two 
invariant lines through the origin, the z-axis and 
the y-axis. 

(d) The y-shear with factor 4 has only one invariant 
line through the origin, the y-axis. 


Solution 1.3 


The matrix ¢ > represents the scaling with 
factors 1 and 3. 


(a) The image of the point (2,0) under this scaling 
is given by 


i. 0\/2\_ f2 
0 3 O}/  \O}?’ 
so the point (2,0) is a fixed point of this scaling. 


(b) Let (0,c) be an arbitrary point on the y-axis. 
The image of the point (0,c) under this scaling 
is given by 


(0 3) (c= (a) 


The point (0, 3c) lies on the y-axis. Since (0,c) is 
an arbitrary point on the y-axis, we have shown 
that every point on the y-axis has its image on 
that axis. Also, as c varies, this image ranges 
over the whole y-axis. Thus the y-axis is an 
invariant line of this scaling. 


ay 


Solution 2.1 
(a) The image under f of the arbitrary point (c,—c) 


on the line y = —2 is given by 
1 2 c\  {(-e 
3 2 —c} c) 
The point (—c,c) is also on the line y = —2, 


and, as c varies, ranges over the whole of that 
line. Thus this is an invariant line of f. 


(b) The images are: 


(3 2)(-1)=(a)=-(4) 
(3 2)(.4)-(3)-4(3), 
(3.5) (4)= (a) 


(c) The image under f of the arbitrary point (c, 3c) 
on the line y = 3x is given by 


(a)U)o 


The point (4c, 6c) is also on the line y = 52, 
and, as c varies, ranges over the whole of that 
line. Thus this is an invariant line of f. 


(d) The images are: 


(3 2) (3)= (a) =4(3) 


Solution 2.2 


Only in part (c) is a full derivation of the 
characteristic equation given. 


(a) The characteristic equation is 
k?+k-2=0. 


This factorises as (k + 2)(k — 1) = 0, so the 
eigenvalues are —2 and 1. 


(b) The characteristic equation is 
k? — 6k + 23 =0. 
Using the quadratic formula, we obtain 


k = (6+ V36 — 92). 


Since 36 — 92 is negative, there are no real 
solutions to the characteristic equation; so this 
matrix has no eigenvalues. 


(c) The characteristic equation is 
k? — (-0.5 + 0.6)k + (—0.3) — (0.26) = 0; 
that is, 
k* — 0.1k — 0.56 = 0. 
Using the quadratic formula, we obtain 
k= 4(0.1+ V2.25) 
_ O.1+1.5 
= 5 
so the eigenvalues are 0.8 and —0.7. 


Solution 2.3 


(a) In Activity 2.2, we saw that the eigenvalues for 
this matrix are —2 and 1, so we use these values 
for k in the eigenvector equation Ax = kx. 


The eigenvector equation with eigenvalue —2, 
-1 2 x x 
(7) G)-G). 
gives the two equations 
—x+2y= —2z, 
v= —2y. 


These equations both reduce to the equation 
x + 2y =0, so we find that the eigenline 
1 


corresponding to the eigenvalue —2 is y = —52. 


Any vector of the form ( i) c#0, is an 
2 
eigenvector for this eigenline, so two such 


: 2 —4 
eigenvectors are “1 and 9) 


The eigenvector equation with eigenvalue 1, 


(1 0) G)=2G) 
1 0 y y)? 
gives the two equations 
—“x+2y=2, 
C= y. 
These equations both reduce to the equation 


y = 4, so the eigenline corresponding to the 
eigenvalue 1 is y = a. 


Any vector of the form (<). c #0, is an 


eigenvector for this eigenline, so two such 


‘ 1 —2 
eigenvectors are 1 and _9 


(b) In Activity 2.2, we saw that the eigenvalues for 
this matrix are 0.8 and —0.7, so we use these 
values for k in the eigenvector equation 
Ax = kx. 


SOLUTIONS TO ACTIVITIES 


The eigenvector equation with eigenvalue 0.8, 
—0.5 1.3 x x 
( 0.2 - (5) =e (5) , 
gives the two equations 


—0.5% + 1.3y = 0.82, 
0.2% + 0.6y = 0.8y. 


These equations both reduce to the equation 
y = 2, so the eigenline corresponding to the 
eigenvalue 0.8 is y = x. 

Any vector of the form (<). c# 0, is an 
eigenvector for this eigenline, so two such 


: 1 3 
eigenvectors are 1 and 3 


The eigenvector equation with eigenvalue —0.7, 
—0.5 1.3 x x 
( 0.2 is) @ ae (5) : 
gives the two equations 


—0.57 + 1.3y = —0.72, 
0.22 + 0.6y = —0.7y. 


These equations both reduce to the equation 
2x + 13y = 0, so the eigenline corresponding to 


the eigenvalue —0.7 is y = — ea. 


Any vector of the form ( : a c#0, is an 
13 
eigenvector for this eigenline, so two such 


: 13 —26 
eigenvectors are 9 and 4): 


Solution 2.4 
The characteristic equation for this matrix is 
k? —3k-4=0. 


This factorises as (k — 4)(k +1) = 0, so the 
eigenvalues are 4 and —1. 


The eigenvector equation with eigenvalue 4, 
—2 3 x x 
(= 5) ()-4G). 
gives the two equations 


—2¢ + 3y = 4a, 


—2x@ + 5y = Ay. 


These equations both reduce to —2x + y = 0, so the 
eigenline corresponding to the eigenvalue 4 is y = 2a. 
An eigenvector for this eigenline is any vector of the 


form ( ) c #0, for example, the vector ( : ), 


The eigenvector equation with eigenvalue —1, 


(3 3) G)-"G) 
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gives the two equations 


—22 + 3y = —2, 


—2@ + 5y = —y. 


These equations both reduce to —x + 3y = 0, so the 
eigenline corresponding to the eigenvalue —1 is 
y= qe. An eigenvector for this eigenline is any 


vector of the form ( ‘). c #0, for example, the 
3 


vector ( : ) : 


Solution 2.5 
(a) The characteristic equation is 
P+k-6=0. 


This factorises as (k + 3)(k — 2) =0, so the 
eigenvalues are —3 and 2. These are the 
diagonal elements of the matrix. 


(b) The characteristic equation for the matrix D is 
k? —(a+dk+ad=0. 


This factorises as (k — a)(k — d) = 0, so the 
eigenvalues are a and d, which are the diagonal 
elements of D. 


Solution 2.6 


The eigenvector equation with eigenvalue —2, 


—2 0 t\__5(% 
0 2 y) y)? 
gives the two equations 


—22 = —2z, 


2y = —2y. 


The first equation holds for every value of xz and the 
second equation holds only for y = 0. So the 
eigenline corresponding to the eigenvalue —2 is y = 0; 
that is, it is the x-axis (as expected). 


The eigenvector equation with eigenvalue 2, 
Co 2) (3)-2G): 

gives the two equations 
—2¢ = 22, 
2y = 2y. 


The second equation holds for every value of y and 
the first equation holds only for z = 0. So the 
eigenline corresponding to the eigenvalue 2 is x = 0; 
that is, it is the y-axis (as expected). 
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Solution 2.7 


(a) The matrix (5 3 has characteristic equation 


k? — 9k = 0, 


which factorises as k(k — 9) = 0. Thus the 
matrix A has eigenvalues 0 and 9. 


(b) The eigenvector equation with eigenvalue 0, 


6 2 t\_g{ 
9 3 y) — \y J? 
gives the two equations 


6x + 2y = 0, 
9x + 38y = 0. 


These equations both reduce to 3x + y = 0, 
which is the eigenline corresponding to the 
eigenvalue 0. This implies that every point on 
the line 3x + y = 0 is mapped to the origin. 


The eigenvector equation with eigenvalue 9, 
6 2 t\_9/* 
9 3) \y y 

gives the two equations 


6x + 2y = 9a, 


9x + 3y = Oy. 


These equations both reduce to 3x — 2y = 0, so 
this is the eigenline corresponding to the 
eigenvalue 9. 


Solution 2.8 
The characteristic equation for this matrix is 
k?4+1=0. 


This equation has no real solutions for k, so this 
matrix has no eigenvalues. 


Solution 2.9 


(a) The matrix A is triangular, so its eigenvalues 
are the elements on the leading diagonal. Thus 
this matrix has only one eigenvalue, k = 1. The 
eigenvector equation with eigenvalue 1, 


CG) 4G), 


gives the two equations x = xz andz+y=y. 
The first equation holds for every value of x 
(and y), but the second equation reduces to 

x = 0. Thus the y-axis is the only eigenline for 
this matrix. 


(b) We saw in Section 1 that a y-shear has only one 
invariant line through the origin, the y-axis, so 
this result agrees with our geometric 
understanding of shears. 


Solution 3.1 


Step 1. The eigenvalues of A are solutions of the 
characteristic equation 


k—k-6=0. 
This equation factorises as (k — 3)(k + 2) =0, so the 
eigenvalues are k = 3 and k = —2. 
3 0 
Step 2. Let D= (; 3) 


Step 3. The eigenvector equation with eigenvalue 3, 


(0) (a) =86s), 


gives the two equations x + 6y = 3a and x = 3y. 
These both reduce to 7 — 3y = 0. Thus y = $x is the 
eigenline for the eigenvalue 3. 


The eigenvector equation with eigenvalue —2 


(0) G)--G): 


gives the two equations 2 + 6y = —2x and x = —2y. 
These both reduce to x + 2y = 0. Thus y = hu is 
the eigenline for the eigenvalue —2. 


o] 


Step 4. For the eigenvalue 3, with eigenline y = $2, 


we can choose the eigenvector a 


For the eigenvalue —2, with eigenline y = —ha, a 


good choice is & ). 


Step 5. Hence we let P = & e 


Step 6. Thus 
1 /-1 —-2 1/1 2 
So eee ac nay 2 
pt=3(7 a)e8 lr 3): 
(We should now have 
i,6) (3 B73’ Wt @ 
1 0/°  \i -tl 0 —2/5\1 -3/?’ 


and this is easily checked, by multiplying out the 
right-hand side.) 


Solution 3.2 


The required matrix powers are: 


sft Qt By 7h BY. 
ce so) 40 iy? 
ec wk. BVP oy: 7 3. Bey. 
AP aa= (3 3 9 10) \39 38)? 


i wee (4 OV/16 0} (64 0 
Di=Dp?=(5 | oo) @ 21" 


SOLUTIONS TO ACTIVITIES 


Solution 3.3 


(9 8) 


II 
a ™~ 

(ee) 
Or 
ran 
ao 
SY 


Thus 
A*=PD‘P“! 
_f(-1 1 81 0\1/-4 1 
= 14 0 16/5 11 
_1/-81 16 -4 1 
5 81 64 11 
_1/f 340 —65 
~ 5 \ 260 145 
7 68 —13 
~ \ 52 29} * 
Solution 3.4 


(a) 


y= 2x 


Figure S.1 


(b) The image P’ is given by 


(3 8) )-()- 


The construction in Figure $.1 agrees with this. 
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Solution 3.5 
—2 3 
—2 5 
—1, with corresponding eigenlines y = 2” and 
1 


(a) (i) By property (a): since 4 > 0, P’ lies on the 
same side of the line y = $x as P; since —1 < 0, 
P’' lies on the opposite side of the line y = 2x 
to P. 


The matrix A = has eigenvalues 4 and 


By property (b): since |4| = 4, the distance from 
P’ to the line y = $x is 4 times the distance 
from P to that line; since | — 1] = 1, the distance 
from P’ to the line y = 2a is equal to the 
distance from P to that line. 


(You should check these descriptions against 
Figure S.1.) 

(ii) The required descriptions are the same as in 
part (i), but with P replaced by P’, and P’ 
replaced by P’’. In particular, P” is on the same 
side of the line y = 42 as P’, but on the 


3 
opposite side of the line y = 2x to P’. 


(b) The image P” is given by 


@Hiones 


Note that (7, 13) is indeed on the opposite side 
of the line y = 2a to P’. 


Solution 4.1 


Let xo = . The second point in the iteration 


2 
3 
sequence is represented by the vector 


m= Am=(3 3) (3)=(i)- 


The third point in the iteration sequence is 
represented by the vector 


x= Am =(5 3) (2) =(is)- 


So the first three points of the iteration sequence are 
(2,3), (8, 12), (32, 48). 


(Note that it is not necessary to calculate these 
points as above. Since the point (2,3) lies on the 
eigenline y = 3x, which has corresponding 
eigenvalue 4, the coordinates of each point in this 
iteration sequence are 4 times the coordinates of the 
previous point.) 
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Solution 4.2 


(a) The initial point (2,3) lies on the eigenline 
y= 3x, with eigenvalue 4. Hence, from the 
table, the iteration sequence moves away from 
(0,0), remaining on the same half of the line 
y = 3a as (2,3). 

(b) The initial point (2,—2) lies on the eigenline 
y = —x, with eigenvalue —1. Hence, from the 
table, the iteration sequence alternates between 
(2,2) and (—2, 2). 


Solution 4.3 


(a) The second, third and fourth points in the 
iteration sequence with initial point (—3,—1) are 
represented by the vectors: 


Cea a. 
Ce eae: 
ESles- es) 


The first four points of the iteration sequence 
are (—3,—1), (—5, —11), (—27, —37) and 
(—101, —155). 


(b) From Example 4.1, the matrix A” is 


i (20 +3(-1)" 2(4)" — a1) | 


5 \.3(4)" —3(-1)”_3(4)" +. 2(-1)” 


Now we calculate x), = A”xo: 


1 ( —8(4)" —7(-1)" 
( —12(4)" + 7(-1)” ) 


To find the 12*" point in this iteration sequence, 
we substitute n = 11 into the expression for x}: 


( =a 1) ) 
Spay 471)" 


re 
x11 = 


—50 331655 


1 
5 
d ee 
5 

_ ( —6 710885 ) 


—10 066 331 


(c) The ratios are as follows (to 6 d.p.): 


Point (z/,, y/,) Ratio yf, /z', 
(3.1) 0.333 333 
(—5,-11) 2.200 000 

(—27, -37) 1.370370 


(—101, —155) 1.534.653 


(—6 710885, —10 066 331) 1.500 001 


It appears that the ratio is tending to 1.5 once 
again. 


Solution 4.4 


By iteration property (a): since 2 > 0, the points of 
the sequence (2p, Yn) lie on the same side of the 
eigenline y = —z as the initial point; since —3 < 0, 
the points alternate between opposite sides of the 
eigenline y = 4a. 


By iteration property (b): since 
max{|2|,| — 3|} =3 > 1, the sequence moves away 
from (0,0). 

Since | — 3| > |2|, the dominant eigenvalue is —3 and 
the dominant eigenline is y = —x. Hence, by iteration 
property (c) (the Dominant Eigenvalue Property) 


Yn 
In 


o] 


lasn oo. 


Thus the long-term behaviour of the sequence may 
be described as follows. The sequence tends in the 
direction of the line y = —2, moving away from the 
origin. It stays on the same side of that line as the 
initial point, the points alternating between opposite 
sides of the line y = 4a. 


The following figure, which you were not expected to 
produce, illustrates the first few terms of the 
sequence. 


T > 

60 
Xa 

(123, —75) 

Se 


Figure S.2 


SOLUTIONS TO ACTIVITIES 


of 


Solutions to Exercises 


Solution 1.1 


(a) 


(b) 


The rotation r; has only one fixed point, the 
origin. Every line through the origin is an 
invariant line. 


The reflection g,/4 has a line of fixed points, the 
line y = x. It has two invariant lines through the 
origin, the line y = x (which is a line of fixed 
points) and the line y = —z. 


The scaling with factors 2 and —3 has only one 
fixed point, the origin. It has two invariant lines 
through the origin, the a-axis and the y-axis. 


The x-shear with factor 1 has a line of fixed 
points, the x-axis. It has only one invariant line 
through the origin, the z-axis (which is a line of 
fixed points). 


Solution 1.2 


(a) 


The matrix which represents the z-shear with 


factor 1 is al To show that every point 


1 
0 1 
on the z-axis is a fixed point, consider an 
arbitrary point on the x-axis, say (c,0). The 


image of the point (c,0) is given by 


()6)-() 


so the point (c,0) is a fixed point. Since (c,0) is 
an arbitrary point on the z-axis, we have shown 
that every point on this axis is a fixed point. 


To show that the line y = x is not an invariant 
line, we show that the image of a point on the 
line y = x is not on the line y = x. Consider the 
point (1,1). The image of the point (1,1) is 
given by 


(01) G) =): 


The point (2,1) does not lie on the line y = x, so 
the line y = x is not an invariant line of this 
shear. (In fact, no point on the line y = x other 
than the origin has its image on the line y = 2.) 
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Solution 2.1 


(a) 


The matrix e 7, has characteristic equation 


k? —8k+7=0. 


This factorises as (k — 7)(k — 1) = 0, so the 
eigenvalues of this matrix are 7 and 1. 


The eigenvector equation with eigenvalue 7, 
6 1 x x 
(5 2) (a) =7G): 
gives the two equations 


62+ y= 72, 


5a + 2y = Ty. 


These equations both reduce to the equation 
y = x, which is the eigenline corresponding to 
the eigenvalue 7. A corresponding eigenvector 


is ; 
1 
The eigenvector equation with eigenvalue 1, 
6 1 x\ («2 
5 2 y) \y)? 
gives the two equations 


62-+y=a, 


5a + 2y = y. 


These equations both reduce to the equation 
y = —52, which is the eigenline corresponding to 
the eigenvalue 1. A corresponding eigenvector 


(3) 


The matrix Ga —0.990 


0 3.396 
matrix, so its eigenvalues are the elements on 
the leading diagonal. Thus this matrix has 
eigenvalues 3.535 and 3.396. 


) is a triangular 


The eigenvector equation with eigenvalue 3.535, 


(" “3)(5)=s99(5) 
gives the two equations 

3.535a — 0.990y = 3.5352, 

3.396y = 3.535y. 
These equations both reduce to the equation 


y = 0, which is the eigenline corresponding to 
the eigenvalue 3.535. A corresponding 


) 
eigenvector is 0 


The eigenvector equation with eigenvalue 3.396, 


(" -2B)(;)-a(5) 
gives the two equations 

3.5352 — 0.990y = 3.3962, 

3.396y = 3.396y. 


The second equation holds for all values of y. 
The first equation reduces to 


0.1392 — 0.990y =0; that is, y= 382, 


which is the eigenline corresponding to the 


eigenvalue 3.396. A corresponding eigenvector 
. ( 990 
iS | 139 J: 


The matrix € i) has characteristic equation 


k? — 5k =0. 


This factorises as k(k — 5) = 0, so the 
eigenvalues of this matrix are 5 and 0. (This 
matrix represents a flattening.) 


The eigenvector equation with eigenvalue 5, 


(a) G)=8G) 
2 4 Yy y)? 
gives the two equations 
x+2y = 52, 
2z + 4y = 5y. 
These equations both reduce to the equation 


y = 2x, which is the eigenline corresponding to 
the eigenvalue 5. A corresponding eigenvector 


re 
ne. 


The eigenvector equation with eigenvalue 0, 


(2 4) ()=0() 
2 4 y y)? 
gives the two equations 
x+2y=0, 
2a + 4y = 0. 


These equations both reduce to the equation 
y= —ta, which is the eigenline corresponding to 
the eigenvalue 0. (Thus all the points on this 
eigenline are mapped to the origin.) A 


: : ? 2 
corresponding eigenvector is 1 


SOLUTIONS TO EXERCISES 


Solution 2.2 


(a) 


The matrix (5 represents a scaling with 


0 3 
factors 2 and 3. It should have two eigenlines, 
the z-axis and the y-axis, with corresponding 
eigenvalues 2 and 3, respectively. 


: : is a diagonal matrix, so its 
eigenvalues are 2 and 3. 


The matrix ( 


The eigenvector equation with eigenvalue 2, 


2 0 r\_5f{% 
0 3 y) -\yy? 
gives the two equations 
2x2 = 22, 
3y = 2y. 
The first equation holds for all values of x, but 
the second equation reduces to y = 0. Thus the 


eigenline corresponding to the eigenvalue 2 is the 
line y = 0, which is the z-axis, as predicted. 


The eigenvector equation with eigenvalue 3, 


2 0 t\_3/% 
0 3 y) —\y)? 
gives the two equations 
2G = 3a; 
by = 3y. 
The second equation holds for all values of y, 
but the first equation reduces to x = 0. Thus the 


eigenline corresponding to the eigenvalue 3 is the 
line x = 0, which is the y-axis, as predicted. 


The matrix ¢ represents an z-shear with 


0 1 
factor 3. It should have one eigenline, the z-axis 
with corresponding eigenvalue 1 (since the 
points on the x-axis are fixed points). 


The matrix is a triangular matrix, so 


1 3 
0 1 
its eigenvalues are on the leading diagonal. Thus 
this matrix has only one eigenvalue, k = 1. 


The eigenvector equation with eigenvalue 1, 


1 3 Z\ [2 
0 1 y) \yJ? 
gives the two equations 
a+ 3y =a, 
yY=y. 
The second equation holds for all values of y, 
but the first equation reduces to y = 0. Thus the 
eigenline corresponding to the eigenvalue 1 is the 
line y = 0, which is the z-axis, as predicted. 


59 


CHAPTER B3 ITERATION WITH MATRICES 


(c) 


. - represents a reflection 
in the line y = —a, which makes an angle of 
3/4 with the positive z-axis. It should have 
two eigenlines. One is the axis of reflection 

y = —« with corresponding eigenvalue 1, since 
this line consists of fixed points. The other is the 
line perpendicular to the axis of reflection, y = x 
with corresponding eigenvalue —1, since points 
on this line are mapped to the same line on the 
opposite side of the origin but remain at the 
same distance from the origin. 


The matrix ( 


The matrix « a) has characteristic 
equation 
k? -1=0. 


This factorises as (k — 1)(k +1) =0, giving the 
two eigenvalues 1 and —1. 


The eigenvector equation with eigenvalue 1, 


(1 0) ()=G): 


gives the two equations 


-y=2, 
—L=y. 
These equations both reduce to y = —z, which is 


the eigenline corresponding to the eigenvalue 1, 
as predicted. 


The eigenvector equation with eigenvalue —1, 


Q: =h t\_ 4/2 
-1 O y) y)’ 
gives the two equations 


-y= a, 

—“£=—y. 
These equations both reduce to y = x, which is 
the eigenline corresponding to the eigenvalue 
—1, as predicted. 
-1 

0 -1 

rotation through an angle of 7 radians about the 
origin. It also corresponds to a uniform scaling 
with factor —1. It should have every line 
through the origin as an eigenline, with 
corresponding eigenvalue —1. 


The matrix ( corresponds to a 


; > is diagonal, so its 
eigenvalues are on the leading diagonal. Thus it 
has only one eigenvalue, k = —1. 


The matrix (~ 


The eigenvector equation with eigenvalue —1, 


Co a) G)=7G)> 


60 


gives the two equations 
-—“£=-2, 
—y=-y. 
Both these equations hold for all values of x and 


y, so every line through the origin is an 
eigenline, as predicted. 


Solution 3.1 
(a) We follow the strategy, using the results from 


Exercise 2.1. 


(i) Step 1. The eigenvalues of this matrix are 7 
and 1. 


7 O 
Step 2. Let D = é i 


Step 3. The eigenlines of this matrix are y = x 
and y = —5a. 


Step 4. An eigenvector for the eigenvalue 7 is 


( : ). An eigenvector for the eigenvalue | is 


1 
5)" 
1 1 
Step 5. Hence we let P = ( ). 


1 —-5 
Step 6. Thus 


1 /-5 -1 1/5 ik 
eh _ ee ds 
Finally, we have 


(5 2)=( 5) (0 1)a(i 1): 


(ii) Step 1. The eigenvalues of this matrix are 
3.535 and 3.396. 


Step 2. Let D= Ca ) 


0 3.396 


Step 3. The eigenlines of this matrix are y = 0 
and y = gana. 


Step 4. An eigenvector for the eigenvalue 3.535 
is ({ . An eigenvector for the eigenvalue 3.396 


. ( 990 
iS | 199 }- 
Step 5. Hence we let P = (; eal 


0 139 
Step 6. Thus 


1 (139 —990 
pola ; 
ama ( 0 1 ) 


Finally, we have 


3.935 —0.990 
0 3.396 


_ (1 990 3.535 0 1 (139 —990 
~ \0 139 0 3.396 / 139 \ 0 it , 


(iii) Step 1. The eigenvalues of this matrix are 5 
and 0. 


5 0 
Step 2. Let D= ( a 


Step 3. The eigenlines of this matrix are y = 2x 
and y = —$a. 


Step 4. An eigenvector for the eigenvalue 5 is 


e} An eigenvector for the eigenvalue 0 is 


2 
aq ye 
1 2 
Step 5. Hence we let P = ( i, 


2 -1 
Step 6. Thus 


1 -1 -2 1/1 2 
ae ae 
Pee aya at) 


Finally, we have 


(2 a)=(2 2) (0 o)s( -1)- 


(b) We multiply the three matrices PDP~' from 


part (a)(i): 
3) (0 alt 4) 


as required. 


(c) First note that 


7? 0 16807 0 
— = 
al eet. Oy 
Then 

A°>=PD°P"! 


_ (1 1) (16807 0\1/5 1 
~ AL =5 0 1/6\1 -1 
_1/16807 1\/5 1 

~ 6\16807 -5/\1 -1 
Ses oo) 


~ 6 \ 84030 16812 


— {14006 2801 
~ \ 14005 2802 ]° 


SOLUTIONS TO EXERCISES 


Solution 3.2 


The matrix A = ( 


: 7 has eigenvalues 7 and 1, 


with corresponding eigenlines y= x and y = —5z. 


(a) 


By property (a): since 7 > 0, P’ lies on the same 
side of the line y = —5a as P; since 1 > 0, P’ 
lies on the same side of the line y = x as P. 


By property (b): since |7| = 7, the distance from 
P’ to the line y = —5z is 7 times the distance 
from P to that line; since |1| = 1, the distance 
from P’ to the line y = z is equal to the distance 
from P to that line. 


The image P’ is given by 
61 1\ (4 
5 2 -2)/ \1/)° 


The locations of P and P’ shown in Figure $8.3 
match the description in part (a). 


Figure 8.3 


As described in part (a): P and P’ lie on the 
same side of the line y = x and on the same side 


of the line y = —5z, they are equidistant from 
y = x, and P’ is 7 times as distant from y = —5a 
as P is. 
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CHAPTER B3 ITERATION WITH MATRICES 


So 


lution 4.1 


1/1 2 
af. 
Bee). 


Th 


en A= PDP“, and 
A” =PD"P"! 


Solution 4.2 


6 1 
5 2 


From Exercise 2.1(a), the matrix A = ( ) has 


eigenvalues 7 and 1 with corresponding eigenlines 
y=aand y= —5a. 


(a) 
(b) 


The point (2,2) lies on the eigenline y = x, and 
the point (1, —5) lies on the eigenline y = —5a. 


The initial point (1,3) does not lie on an 
eigenline of A. So in part (i), the iteration 
properties of generalised scalings are used. The 
other initial points lie on the eigenlines of A. So 
in parts (ii) and (iii) the table on page 41 is used. 


(i) By iteration property (a): since 7 > 0, the 
points of the sequence (2n, Yn) with initial point 
(1,3) lie on the same side of the eigenline 

y = —5z as the initial point; since 1 > 0, the 
points lie on the same side of the eigenline y = z. 


By iteration property (b): since 
max{|7|,|1|} = 7 > 1, the sequence moves away 
from (0,0). 


Since |7| > |1], the dominant eigenvalue is 7 and 
the dominant eigenline is y= «. Hence by 
iteration property (c) (the Dominant Eigenvalue 
Property) 

am og lasn-o. 

In 
Thus the sequence tends in the direction of the 
line y = x, moving away from the origin. It stays 
on the same side of that line as the initial point 
and also on the same side of the line y = —5z. 


(ii) The initial point (2,2) lies on the eigenline 
y = , with eigenvalue 7. Hence, from the table, 
the iteration sequence moves away from (0,0), 
remaining on the same half of the line as (2, 2). 
(iii) The initial point (1,—5) lies on the 
eigenline y = —5a, with eigenvalue 1. Hence, 
from the table, the iteration sequence is the 


constant sequence (2p, Yn) = (1, —5), for 
i ae 
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